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Abstract 

Recently, several authors have studied maps where a function, describing the local diffu- 
sion matrix of a diffusion process with a linear drift towards an attraction point, is mapped 
into the average of that function with respect to the unique invariant measure of the dif- 
fusion process, as a function of the attraction point. Such mappings arise in the analysis 
of infinite systems of diffusions indexed by the hierarchical group, with a linear attractive 
interaction between the components. In this context, the mappings are called renormaliza- 
tion transformations. We consider such maps for catalytic Wright-Fisher diffusions. These 
are diffusions on the unit square where the first component (the catalyst) performs an au- 
tonomous Wright-Fisher diffusion, while the second component (the reactant) performs a 
Wright-Fisher diffusion with a rate depending on the first component through a catalyzing 
function. We determine the limit of rescaled iterates of renormalization transformations 
acting on the diffusion matrices of such catalytic Wright-Fisher diffusions. 
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Part I 

1 Introduction and main result 

Several authors |BCGdH95! IBCGdH971 ldHS981 ISch981 ICDG04| have studied maps where 
a function, describing the local diffusion matrix of a diffusion process, is mapped into the 
average of that function with respect to the unique invariant measure of the diffusion process 
itself. Such mappings arise in the analysis of infinite systems of diffusion processes indexed by 
the hierarchical group, with a linear attractive interaction between the components D G93al 
DG96, DGV95]. In this context, the mappings are called renormalization transformations. We 
follow this terminology. For more on the relation between hierarchically interacting diffusions 
and renormalization transformations, see Appendix I A. I! 

Formally, such renormalization transformations can be defined as follows. 

Definition 1.1 (Renormalization class and transformation) Let D C M. d be nonempty, 
convex, and open. Let W be a collection of continuous functions w from the closure D into the 
space Mi of symmetric non-negative definite d x d real matrices, such that Xw £ W for every 
A > 0, w £ W. We call W a prerenormalization class on D if the following three conditions 
are satisfied: 

(i) For each constant c > 0, w £ W, and x £ D, the martingale problem for the operator 
Ax™ is well-posed, where 

d d 

Arm ■■= E c (^ - w)&f(v) + E ^M^m (y e d), (i.i) 

i=l i,j=l 

and the domain of Ax' w is the space of real functions on D that can be extended to a 
twice continuously differentiable function on M rf with compact support. 

(ii) For each c > 0, w £ W, and x £ D, the martingale problem for A x ' has a unique 
stationary solution with invariant law denoted by Vx ' ■ 

(iii) For each c > 0, w £ W, x £ D, and i,j = l,...,d, one has / (dy)\wij (y)\ < oo. 

J~D 

If W is a prerenormalization class, then we define for each c > and w £ W a matrix- valued 
function F c w on D by 



F c w(x) ■= _"x w (dy)w(y) (x e D). (1.2) 
Jd 

We say that W is a renormalization class on D if in addition: 
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(iv) For each c > and w G W, the function F c w is an element of W. 

If W is a renormalization class and c > 0, then the map F c : W — > W defined by (|1.2j) 
is called the renormalization transformation on W with migration constant c. In p. 1)1 . io is 
called the diffusion matrix and x the attraction point. O 

Remark 1.2 (Associated SDE) It is well-known that -D-valued (weak) solutions y = 
(y 1 , . . . ,y d ) to the stochastic differential equation (SDE) 

n 

dyi = c(^-yj)dt + V2^ay(y t )d^ (t > 0, i = 1, . . . ,d), (1.3) 

i=i 

where -B = (B l ,...,B n ) is n-dimensional (standard) Brownian motion (n > 1), solve the 
martingale problem for A x ' w if the d x n matrix-valued function a is continuous and satisfies 
J2k °ik°~3k = Wij. Conversely EK86, Theorem 5.3.3], every solution to the martingale problem 
for A x ' w can be represented as a solution to the SDE 1)1.3)1 . where there is some freedom in 
the choice of the root a of the diffusion matrix w. O 

In the present paper, we concern ourselves with the following renormalization class on [0, l] 2 . 



Definition 1.3 (Renormalization class of catalytic Wright-Fisher diffusions) We set 

>V cat : = { W ^P : a > 0, p G Tt}, where 

^):=( aXl{1 - Xl) p(ll)l2 ° (1 _ X2) ) (0,1]'), (1.4) 

and 

TL := {p : p a real function on [0, 1], p > 0, p Lipschitz continuous}. (1-5) 
Moreover, we put 

7-Ll,r :={peH: l{ p (o)>o} = I, 1 fr(l)>0} = r ) r = 0, 1), (1.6) 

and set W^t := {w a ' p : a > 0, p £ Hi >r } (I, r = 0, 1). O 

By B,emark ll.21 solutions y = (y 1 , y 2 ) to the martingale problem for A x ' w can be represented 
as solutions to the SDE 

(i) dy, 1 = c(xi - y^dt + J2ayj(l - y])dB t \ 

V ; (1-7) 

(h) dy 2 = c (x 2 - y 2 )dt + ^(yj)yf(l - y 2 )di? 2 . 

We call y 1 the Wright-Fisher catalyst with resampling rate a and y 2 the Wright-Fisher reactant 
with catalyzing function p. 

For any renormalization class W and any sequence of (strictly) positive migration constants 
(cfc)fc>o, we define iterated renormalization transformations F^> : W — ► W, as follows: 

F( n+1 )u; := F Cn (F^w) (n > 0) with F<®w:=w (w G W C at)- (1.8) 
We set so := and 

n— 1 - 

s n :=V— (l<n<oo). (1.9) 

^— ' Cfc 

fc=0 ft 

Here is our main result: 
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Theorem 1.4 (Main result) 

(a) The set W ca t is a renormalization class on [0, l] 2 and F c (W c £ t ) C W c £ t (c > 0, l,r = 0, 1). 

(b) Fix (positive) migration constants (ck)k>o such that 

(i) s n — > oo and (ii) n+1 — > 1 + 7* (1-10) 

n— »oo s n n— >oo 

for some 7* > 0. If w € Wj£± {1,t = 0, 1), i/ien uniformly on [0, l] 2 , 

s n F( n )w — ► «;*, (1.11) 

n^oo 

where the limit w* is the unique solution in to the equation 

(i) (1 + 7*)^*™* = ™* 

2 

(") \Y1 W U X )S^ W ^ X ) +W ^ X )= {] (^[0,1] 2 ) 

(c) The matrix w* is of the form w* = w l,p , where p* = p\ ri * S Tli :r depends on I, r, and 7*. 
One has 

P0A7* = and 7 » = 1 /or a// 7* > 0. (1-13) 

For eac/i 7* > 0, the function p^ 1 » is concave, nondecreasing, and satisfies p$ 1 » (0) = 0, 
Poi 7 *(l) = By symmetry, analoguous statements hold for p\ 7 * ■ 

Conditions 1)1. 10|) (i) and (ii) are satisfied, for example, for Ck = (1 +7*) -fc . Note that the 
functions ^0,0,7* an d Pi, 1,7* are independent of 7* > 0. We believe that on the other hand, 
Pq 1 » is not constant as a function of 7*, but we have not proved this. 
The function 1 i s the unique nonnegative solution to the equation 

l x (l- x )£ I p( x ) + p( x )(l-p( x )) = (x€[0,l]) (1.14) 

with boundary conditions p(0) = and p(l) > 0. This function occurred before in the work of 
Greven, Klenke, and Wakolbinger GKVY01 formulas (1.10)— (1.11)]. In Section l4~Tl we discuss 
the relation between their work and ours. 

Outline In Part I of the paper (Sections we present our results and our main techniques 
for proving them. Part II (Sections contains detailed proofs. Since the motivation for 
studying renormalization classes comes from the study of linearly interacting diffusions on the 
hierarchical group, we explain this connection in Appendix lAl 

Outline of Part I In the next section, we place our main result in a broader context. We 
give a more thorough introduction to the theory of renormalization classes on compact sets 
and discuss earlier results on this topic. In Section EJ we discuss special properties of the 
renormalization class W ca t from Definition II. HI In particular, we show how techniques from 
the theory of spatial branching processes can be used to prove Theorem 11.41 In Section 21 we 
discuss the relation of our work with that in |GKW0i] and mention some open problems. 

Notation If E is a separable, locally compact, metrizable space, then C(E) denotes the space 
of continuous real functions on E. If E is compact then we equip C(F) with the supremumnorm 



if 7* > 0, 



if 7* = 0. 



(1.12) 
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II • I |oo- We let B(E) denote the space of all bounded Borel measurable real functions on E. 
We write C+(E) and C [0A] (E) for the spaces of all / G C(E) with / > and < / < 1, 
respectively, and define B + (E) and B[ ^(E) analogously. We let Ai(E) denote the space of 
all finite measures on E, equipped with the topology of weak convergence. The subspaces of 
probability measures is denoted by A4i(E). We write M(E) for the space of finite counting 
measures, i.e., measures of the form v = YliLi with x±, . . . , x m G E (m > 0). We interpret 
v as a collection of particles, situated at positions x±, . . . ,x m . For fieM(E) and / G B(E) 
we use the notation (n, f) := J E f d/z and |/i| := fJ-(E). By definition, T>e[0,oo) is the space 
of cadlag functions w : [0, oo) — » E, equipped with the Skorohod topology. We denote the law 
of a random variable y by C{y). If y = (yt)t>o is a Markov process in E and x G E, then 
P x denotes the law of y started in yo = x. If [i is a probability law on E then P^ denotes 
the law of y started with initial law £(yo) = [i- For time-inhomogeneous processes, we use 
the notation P l ' x or P to denote the law of the process started at time t with initial state 
yt = x or initial law C{yt) = M> respectively. We let E x ,E fl , . . . etc. denote expectation with 
respect to P X ,P IJ ', . . ., respectively. 

2 Renormalization classes on compact sets 
2.1 Some general facts and heuristics 

In this section, we explain that our main result is a special case of a type of theorem that 
we believe holds for many more renormalization classes on compact sets in Moreover, we 
describe some elementary properties that hold generally for such renormalization classes. The 
proofs of Lemmas 12. 1H2.8I can be found in Section f5.ll below. 

Fix a prerenormalization class W on a set D where D C M d is open, bounded, and convex. 
Then W is a subset of the cone C{D,M+) of continuous M^-valued functions on D. We 
equip C(D,Mf) with the topology of uniform convergence. Our first lemma says that the 
equilibrium measures Ux' w and the renormalized diffusion matrices F c w(x) are continuous in 
their parameters. 

Lemma 2.1 (Continuity in parameters) 

(a) The map (x,c,w) i— > Vx W from D x (0, oo) x W into M.\{D) is continuous. 

(b) The map (x,c,w) i— > F c w(x) from D x (0, oo) x W into Mf is continuous. 

In particular, x i— > Vx W is a continuous probability kernel on D, and F c w G C(D,M^) for all 
c > and w G W. Recall from Definition [TT] that Xw G W for all w G W and A > 0. The 
reason why we have included this assumption is that it is convenient to have the next scaling 
lemma around, which is a consequence of time scaling. 

Lemma 2.2 (Scaling property of renormalization transformations) One has 





(2.1) 



The following simple lemma will play a crucial role in what follows. 
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□ □ □ 



w G W, 



1.1 



W G W, 



0,1 



M)6W 



1.0 



W G W,' 



0,0 



cat ^ r K cat ^ r K cat ^_ >• x cat 

Figure 1: Effective boundaries for w G W ca t. 



Lemma 2.3 (Mean and covariance matrix) For all x £ D and i,j = 1, . . . , d, the mean 
and covariances of v x w are given by 



^ w (dy)(y l -x l )=0, 

D 

(ii) _ u x W ( d y)(yi - x i)(Vj ~ x j) = \F c Wij{x). 
Jd 

For any w £C(D,Mf), we call 

d w D := {x G D : Wij(x) = j = 1, . . . , d} 



(2.2) 



(2.3) 



the effective boundary of D (associated with w). If y is a solution to the martingale problem 
for the operator Yli j=i w ij(v) dy^dy (i' e '' operator in without the drift), then, by 

martingale convergence, y^ converges a.s. to a limit Voo! it is not hard to see that yoo G d w D 
a.s. The next lemma says that the effective boundary is invariant under renormalization. 



Lemma 2.4 (Invariance of effective boundary) One has 8f cW D = d w D for all w G W, 

c> 0. 



For example, for diffusion matrices w from the renormalization class W = W ca t> there 
occur four different effective boundaries, depending on whether w G W^ t , W° a J, W^? , or >V C ^. 
These effective boundaries are depicted in Figure^ The statement from Theorem II .41 (a) that 
F c (w'a t ) C W' at 1S J us t the translation of Lemma |2~11 to the special set-up there. 

From now on, let W be a renormalization class, i.e., W satisfies also condition (iv) from 
Definition 11.11 Fix a sequence of (positive) migration constants (ck)k>o- By definition, the 
iterated probability kernels K w, ^ n "> associated with a diffusion matrix w G W (and the constants 
( c fc)fc>o) are the probability kernels on D defined inductively by 

^' (n+1) (dz):=jT_^' F( " ) -(dy)^W(dz) (n>0) with K^°\dy) := 8 x (dy), (2.4) 
with F( n ) as in (fTHj) . Note that 

F^w{x) = l_K^ n \dy)w{y) (x G D, n> 0). (2.5) 
Jd 

The next lemma follows by iteration from Lemmas 12 . 1 1 and 12 . 31 It their essence, this lemma 
and Lemma 12.61 below go back to BCGdH95 . 
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Lemma 2.5 (Basic properties of iterated probability kernels) For each w £ W, the 
j(w,(n) are con ti nuous probability kernels on D. Moreover, for all x & D, i, j = 1, . . . , d, and 
n > 0, the mean and covariance matrix of K™'^ are given by 

(i) [_K^ n \dy)( yi - Xi )=0, 

Jd (2.6) 
K^ n \dy)( yi - Xi ){y 3 - xj) = s n F^ Wij (x). 



D 

We equip the space C{D,M.\{D)) of continuous probability kernels on D with the topology 
of uniform convergence (since Aii(D) is compact, there is a unique uniform structure on 
.Mi(.D) generating the topology). For 'nice' renormalization classes, it seems reasonable to 
conjecture that the kernels K w >( n ) converge as n —* oo to some limit K w '* in C(D, M.\{D)). If 
this happens, then formula 1)2.6)) (ii) tells us that the rescaled renormalized diffusion matrices 

F (n) 

w converge uniformly on D to the covariance matrix of K w '*. This gives a heuristic 
explanation why we need to rescale the iterates F^ n 'w with the scaling constants s n from 
to get a nontrivial limit in (jl.llj) . 

We now explain the relevance of the conditions (jl.K)j) (i) and (ii) in the present more 
general context. If the iterated kernels converge to a limit K w '*, then condition (jl.lUJI (i) 
guarantees that this limit is concentrated on the effective boundary: 

Lemma 2.6 (Concentration on the effective boundary) If s n — > oo, then for any 

n— >oo 

/ G C(D) such that f = on d w D: 

hm sup / K^ n \dy)f(y) =0. (2.7) 

In fact, the condition s n — > oo guarantees that the corresponding system of hierarchically 
interacting diffusions with migration constants (ck)k>o clusters in the local mean field limit, 
see DG93a, Theorem 3] or Add endix I A . 1 1 b elow . 

To explain also the relevance of condition 1)1.10(1 (ii), we observe that using Lemma l2.2( we 
can convert the rescaled iterates s n F^ n ' into (usual, not rescaled) iterates of another transfor- 
mation. For this purpose, it will be convenient to modify the definition of our scaling constants 
s n a little bit. Fix some j3 > and put 

s n :=(3 + s n (n>0). (2.8) 

Define rescaled renormalization transformations F^ : W — » W by 

F 7 u> := (1 + -y)F 1/7 w (7 > 0, we W). (2.9) 

Using ()2.1(l (ii), one easily deduces that 

s n F in) w = F 7n _ x o • • • oF J0 (Pw) {weW, n> 1), (2.10) 

where 

7n := J— in > 0). (2.11) 

Stj ,C n 
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We can reformulate the conditions (|1,1U|) (i) and (ii) in terms of the constants (j n )n>0- 
Indeed, it is not hard to check 1 that equivalent formulations of condition (jl.lOj) (i) are: 



s n — > OO, (11) s n — > 

n— too. n — >oo 



oo, (iii) ^ 7n = oo. (2.12) 



Since s n+ \/ s n = 1 + 7 n we see moreover that, for any 7* £ [0, 00], equivalent formulations of 
condition (|1.1U|) (ii) are: 

(i) f2±l _ l + 7 *, (ii) l +7 *, (iii ) 7n _^ 7 *. ( 2 .13) 

S n n— >oo s n n—foo n—*co 

If < 7* < 00, then, in the light of l|2.10j) . we expect s n F^w to converge to a fixed point 
of the transformation F 7 *. If 7* = 0, the situation is more complex. In this case, we expect 
the orbit s n F^ n 'w 1— * s n+ iF^ n+1 ^w 1— * ■ ■ ■ , for large n, to approximate a continuous flow, the 
generator of which is 

d 

^ 1 -^(F 1 w-wYx) = \Y^w ij {x)Q^-w{x) + w{x) (x G ~D). (2.14) 



To see that the right-hand side of this equation equals the left-hand side if w is twice contin- 
uously differentiable, one needs a Taylor expansion of w together with the moment formulas 
(|2.2|) for Ux^" l ' w . Under condition condition (|2.12l) (iii), we expect this continuous flow to reach 
equilibrium. 

In the light if these considerations, we are led to at the following general conjecture. 

Conjecture 2.7 (Limits of rescaled renormalized diffusion matrices) Assume that 
s n — > 00 and s n+ i/s n — > 1 + 7* for some 7* G [0, 00]. Then, for any w G W, 



s n F^w — ► w*, (2.15) 

n— >oo 

where w* satisfies 

(i) T^*w* =w* if < 7* < 00, 

d 

(ii) l 2 Y,w* J (x)^-w*(x)+w*(x)=0 (xeD) ? / 7 * = 0, (2.16) 

(iii) lim F^w* = w* if 7* = 00. 

7— »oo 

We call (|2.1(i|) (ii), which is in some sense the 7* — > limit of the fixed point equation 
(|2.16j) (i), the asymptotic fixed point equation. A version of formula ()2.16|) (ii) occured in 
Swa99, formula (1.3.5)] (a minus sign is missing there). 

In particular, one may hope that for a given effective boundary, the equations in (|2.16|) 
have a unique solution. Our main result (Theorem I1.4JI confirms this conjecture for the 
renormalization class W ca t and for 7* < 00. In the next section, we discuss numerical evidence 
that supports Conjecture 12.71 in the case 7* = for other renormalization classes on compacta 
as well. 



x To see this, let Soo G (0, 00] denote the limit of the s„ and note that on the one hand, ^2 n l/(s n c„) > 
^„log(l + l/(s„c n )) = log(IT n s n+ i/s„) = loglXo/si), while on the other hand £ n l/(s n c„) < n„( 1 + 

1/(S„C„)) =Il„Sn + l/Sn =Soo/si. 



s 



In previous work on renormalization classes, fixed shapes have played an important role. 
By definition, for any prerenormalization class W, a fixed shape is a subclass W C W of the 
form W = {Xw : A > 0} with / w £ W, such that F C (VV) C W for all c > 0. The next 
lemma describes how fixed shapes for renormalization classes on compact sets typically arise. 

Lemma 2.8 (Fixed shapes) Assume that for each < 7* < 00, there is a 7^ w* = u>*» £ 
W such that s n F^w — > w** whenever w 6 W, s n — » oo ; and s n+ i/s n — ► 1 + 7*. Then: 

(a) «;** is i/ie unique solution in W 0/ equation \2.16\) (i). 

(b) //w* = does no£ depend on 7*, i/ien 

Fc(Aw * ) = ( i +I)-V (A, c > 0). (2.17) 
Moreover, {Xw* : A > 0} is the unique fixed shape in W. 

(c) If the w** for different values of 7* are not constant multiples of each other, then W 
contains no fixed shapes. 

Note that by Theorem O is a renormalization class satisfying the general assumptions 
of Lemma 12.81 The unique solution of ()2.16j) (i) in W^t is of the form w* = w 1,p * where 
p* = Pq 1 » . We conjecture that the Pq t » for different values of 7* are not constant multiples 

of each other, and, as a consequence, that W^t contains no fixed shapes. 

Many facts and conjectures that we have discussed can be generalized to renormalization 
classes on unbounded D, but in this case, the second moments of the iterated kernels K w '^ 
may diverge as n — > 00. As a result, because of formula (|2,6|) (ii), the s n may no longer be 
the right scaling factors to find a nontrivial limit of the renormalized diffusion matrices; see, 
for example, |B(X4dH97j . 

2.2 Numerical solutions to the asymptotic fixed point equation 

Let 1 1 — > w(t, • ) be a solution to the continuous flow with the generator in (|2.14|) . i.e., w is an 
M^-valued solution to the nonlinear partial differential equation 

d 

§tw(t,x) = \ mj(t,x)g^:w(t,x) + w(t,x) (t>0,xeD). (2.18) 

Solutions to (|2.18|) are quite easy to simulate on a computer. We have simulated solutions 
for all kind of diffusion matrices (including nondiagonal ones) on the unit square [0, l] 2 , with 
the effective boundaries 1-6 depicted in Figure |2j For all initial diffusion matrices w(0, ■ ) we 
tried, the solution converged as t — > 00 to a fixed point w*. In all cases except case 6, the 
fixed point was unique. The fixed points are listed in Figure EJ The functions 1 an< ^ 1* 
from Figure |^1 are plotted in Figure |31 Here p$ 1 is the function from Theorem 11.41 (c) . 

The fixed points for the effective boundaries in cases 1,2, and 4 are the unique solutions 
of equation (|1.12j) (ii) from Theorem 11.41 in the classes W^, W^ti an d VV^t' respectively. 
The simulations suggest that the domain of attraction of these fixed points (within the class 
of "all" diffusion matrices on [0, l] 2 ) is actually a lot larger than the classes W c ^ t , W c ^ t , and 
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case 


effective boundary 


fixed points w* of (j2.18j) 


1 
2 
3 
4 
5 
6 




/xi(l-xi) \ 
V x 2 (l-x 2 )J 

V Po,l,o(xi)x 2 (l-X2)) 

(q*(xi,x 2 ) o \ 

V q*{x 2 , Xl )J 
f Xl (l- Xl ) 0\ 

V o V 

- xi)l {x2> o} 0\ 

g{xi,x 2 )[ 

\m 21 m 22 J 



Figure 2: Fixed points of the flow ()2.18jl . 



The function q* from case 3 satisfies q*{x\, 1) = x\{l — x\) and is zero on the other parts 
of the boundary. In contrast to what one might perhaps guess in view of case 2, q* is not of 
the form q*(xi,x 2 ) = f{x 2 )xi(l — x\) for some function /. 

Case 5 is somewhat degenerate since in this case the fixed point is not continuous. 

The only case where the fixed point is not unique is case 6. Here, m can be any positive 
definite matrix, while g*, depending on m, is the unique solution on (0, l) 2 of the equation 
1 + | Ylfj=i m ij dx dx 9*( x ) = 0> with zero boundary conditions. 

2.3 Previous rigorous results 

In this section we discuss some results that have been derived previously for renormalization 
classes on compact sets. 

Theorem 2.9 |BCGdH95l IDGV95 | (Universality class of Wright-Fisher models) 

Let D := {x £ M. d : X{ > Vi, Y^i=i x i < 1}> an< ^ ^ { e 0> • • • > e d}, with e$ := (0, . . . , 0) and 
ei := (1, 0, . . . , 0), . . . , ed := (0, . . . , 0, 1) be the extremal points of D. Let w*^{x) := Xi(5ij — xj) 
(x G D, i, j = 1, . . . , d) denote the standard Wright-Fisher diffusion matrix, and assume that 
W is a renormalization class on D such that w* G W and d w D = {eo, • • • , e^} for all w G W. 
Let (ck)k>o be migration constants such that s n — > oo as n — > oo. Then, for all w G W, 
uniformly on D, 

s n F^w — ► w*. (2.19) 

The convergence in (|2.19|) is a consequence of Lemmas 12.51 and 12.61 The first moment formula 
(US!! (i) and (J22J) show that K^' (n> converges to the unique distribution on {eo, • • • , e^} with 
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Figure 3: The functions Pq 1 o anc ^ °* f rom cases 2 and 3 of Figured 

mean x, and by the second moment formula Q2.6|) (ii) this implies the convergence of s n F^w. 

In order for the iterates in (|2,19|) to be well-defined, Theorem 12. 91 assumes that a renormal- 
ization class W of diffusion matrices w on D with effective boundary {eo, • • • , e^} is given. The 
problem of finding a nontrivial example of such a renormalization class is open in dimensions 
greater than one. In the one-dimensional case, however, the following result is known. 

Lemma 2.10 | DG9 3bj (Renormalization class on the unit interval) The set 

W DG := {w G C[0, 1] : w = on {0, 1}, w > on (0, 1), w Lipschitz} (2.20) 
is a renormalization class on [0, 1]. 

About renormalization of isotropic diffusions, the following result is known. Below, dD := 
D\D denotes the topological boundary of D. 

Theorem 2.11 |dHS98j (Universality class of isotropic models) Let D C M. d be open, 
bounded, and convex and let m G be fixed and (strictly) positive definite. Set w*j(x) := 

m ij9*( x )j where g* is the unique solution of 1 + 2 m ij dx dx 9* ( x ) = ^ f or x ^ D an< ^ 
g*(x) = for x £ dD. Assume that W is a renormalization class on D such that w* G W 
and such that each w £ W is of the form 

Wij(x) = rriijg(x) (x £ D, i,j = l,...,d), (2-21) 

for some g G C(D) satisfying g > on D and g = on dD. Let (ck)k>o be migration constants 
such that s n — > 00 as n — > 00. Then, for all w G W, uniformly on D, 

s n F {n) w — ► w*. (2.22) 

n— >oo 

The proof of Theorem 12.111 follows the same lines as the proof of Theorem 12 .9( with the 
difference that in this case one needs to generalize the first moment formula (|2.6|) (i) in the 
sense that J-pK^idyWy) = h(x) for any m-harmonic function h, i.e., h G C(D) satisfying 
mi i dx^dx h{ x ) = f° r x £ D. The kernel K™'^ now converges to the m-harmonic 
measure on dD with mean x, and this implies (|2.22|) . 

Again, in dimensions d > 2, the problem of finding a 'reasonable' class W satisfying the 
assumptions of Theorem 12. Ill is so far unresolved. The problem with verifying conditions (i)- 
(iv) from Definition 1 1.1 1 in an explicit set-up is that (i) and (ii) usually require some smoothness 
of w, while (iv) requires that one can prove the same smoothness for F c w, which is difficult. 
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The proofs of Theorems 12.91 and 12.111 are based on the same principle. For any diffusion 
matrix w, let H w denote the class of w-harmonic functions, i.e., functions h E C(D) satisfying 
Tlij Wi 3 ( x ) dx dx M x ) = on D. If w belongs to one of the renormalization classes in Theo- 
rems EH1 and then H w has the property that T£ t h(H w ) C H w for all c > 0, x £ D, and 

t > 0, where T^ t h(y) := h(x + (y — x)e~ ct ) is the semigroup with generator Yli=i c ( x i ~Vi)^p^ 
i.e., the operator in (|l.ljl without the diffusion part. In this case we say that w has invariant 
harmonics; see SwaOO . As a consequence, one can prove that the iterated kernels satisfy 
fjy Kx^ n \dy)h{y) = h(x) for all h € H w and If s n — > oo, then this implies that K™ 

converges to the unique ff^-harmonic measure on d w D with mean x. Diffusion matrices from 
W cat do not in general have invariant harmonics. Therefore, to prove Theorem 11.41 we need 
new techniques. 

Note that in the renormalization classes from Theorems l2.9l and l2.lll the unique attraction 
point w* does not depend on 7*. Therefore, by Lemma 12.81 these renormalization classes 
contain a unique fixed shape, which is given by {Xw* : A > 0}. 

3 Connection with branching theory 

From now on, we focuss on the renormalization class W ca t- We will show that for this renor- 
malization class, the rescaled renormalization transformations Fj from (|2.9|) can be expressed 
in terms of the log-Laplace operators of a discrete time branching process on [0, 1]. This will al- 
low us to use techniques from the theory of spatial branching processes to verify Coniecture l2.7l 
for the renormalization class W ca t in the case 7* < 00. 

3.1 Poisson-cluster branching processes 

We first need some concepts and facts from branching theory. Finite measure- valued branching 
processes (on IR) in discrete time have been introduced by Jifina |.Tir64j . We need to consider 
only a special class. Let E be a separable, locally compact, and metrizable space. We call a 
continuous map Q from E into Mi{M{E)) a continuous cluster mechanism. By definition, 
an M (E)- valued random variable X is a Poisson cluster measure on E with locally finite 
intensity measure [i and continuous cluster mechanism Q, if its log-Laplace transform satisfies 

-log£[e-<*'/>] = //x(dx)(l- / Q(x,dx)e-<X,/>) (/ 6 B + (E)). (3.1) 
Je v Jm{e) / 

For given /i and Q, such a Poisson cluster measure exists, and is unique in distribution, 
provided that the right-hand side of <|3.1|) is finite for / = 1. It may be constructed as X = 
J2i Xxn where ^ 5 Xi is a (possibly infinite) Poisson point measure with intensity fi, and given 
xi, X2, ■ ■ ., the Xxn Xx 2 i ■ ■ ■ are independent random variables with laws Q(xi, • ), Q(x2, •)>•••> 
respectively. 

Now fix a finite sequence of functions £ C + (E) and continuous cluster mechanisms Qk 
(k = 1, . . . , n), define 

U k f(x):=q k (x)(l- [ Q h (x,d X )e-<X>f)) (x € E, f G B+{E), k = 1, . . . , n), (3.2) 
v Jm(e) j 

and assume that 

sup^l^) < 00 [k = 1, . . . , n). (3-3) 

x&E 
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ThenUk maps B + {E) into B + (E) for each k, and for each A^(£')-valued initial state Xq, there 
exists a (time-inhomogeneous) Markov chain (Xq, . . . , X n ) in A4(E), such that X^, given Xk-i, 
is a Poisson cluster measure with intensity q^X^-i and cluster mechanism Qk- It is not hard 
to see that 

£^ e -(* n ,/)j = e -(/i,Wio---oW„/) (fj,eM(E), f e B + (E)). (3.4) 

We call <Y = (A?o, . . . , <-t n ) the Poisson- cluster branching process on E 1 with weight functions 
qi, . . . ,q n and cluster mechanisms Q 1} . . . , Q n . The operator is called the log-Laplace op- 
erator of the transition law from A^-i to Xk- Note that we can write (|3.4j) in the suggestive 
form 

[Pois(/* n ) = 0] = P[Pois((«i o • • • o U n f)v) = 0] . (3.5) 

Here, if /i is an M (P)-valued random variable, then Pois(^i) denotes an M (E)-valued random 
variable such that conditioned on /x, Pois(/x) is a Poisson point measure with intensity /j,. 

3.2 The renormalization branching process 

We will now construct a Poisson-cluster branching process on [0, 1] of a special kind, and show 
that the rescaled renormalization transformations on W ca t can be expressed in terms of the 
log-Laplace operators of this branching process. 

By Lemma 15.41 below, for each 7 > and x G [0, 1], the SDE 

dy(t) = i {x - y(t))dt + >/2y(t)(l - y(t))dB(t), (3.6) 

has a unique (in law) stationary solution. We denote this solution by (y2(£))teR- Let r 7 be 
an independent exponentially distributed random variable with mean 7, and set 

Z2 := r S y2( _ t/2) dt (7 > 0, x G [0, 1]). (3.7) 

JO 

Define constants g 7 and continuous (by Corollary I5.1UI below) cluster mechanisms Q 7 by 

? 7 :=i + l and Q 7 (x, • ) := (7 > 0, x G [0, 1]), (3.8) 

and let W 7 denote the log-Laplace operator with (constant) weight function g 7 and cluster 
mechanism Q 7 , i.e., 

Uyf(x):=qJl- f Q 7 (x,d X )e-(xJ}) (x G [0, 1], / G P+[0, 1], 7 > 0). (3.9) 
v JM([o,i\) ' 

We now establish the connection between renormalization transformations on W ca t and log- 
Laplace operators. 

Proposition 3.1 (Identification of the renormalization transformation) Let F 7 be 

the rescaled renormalization transformation on W ca t defined in i2.9\l . Then 

F^w 1 ^ = w 1 > U ~f p {pen, 7>0). (3.10) 
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Fix a diffusion matrix w a ' p 6 W ca t and migration constants (ck)k>o- Define constants s n and 
7 n as in ()2.8)) and 1)2.11)) . respectively, where /3 := 1/a. Then Proposition 13.11 and formula 
()2~Tn|l show that 

s nJ F( n )w Q ' p = wWn-i °---°W 70 (f)_ (3 1X ) 

Here W 7n _ 1 , . . . ,W 7o are the log-Laplace operators of the Poisson-cluster branching process X = 
{X- n , . . . , Xq) with weight functions q ln _ 1 , . . . ,q lQ and cluster mechanisms Q 7n _ i; • • • , Qj - 
We call X (started at some time — n in an initial law C(X- n )) the renormalization branching 
process. By formulas ()3.4|) and 1)3.11)1 . the study of the limiting behavior of rescaled iterated 
renormalization transformations on W C at reduces to the study of the renormalization branching 
process X in the limit n —* oo. 



3.3 Convergence to a time- homogeneous process 

Let X = (X_ n , . . . , Xq) be the renormalization branching process introduced in the last section. 
If the constants (7fc)fc>0 satisfy ^n7« = oo and 7„ — > 7* for some 7* £ [0, 00), then X is 
almost time-homogeneous for large n. More precisely, we will prove the following convergence 
result. 

Theorem 3.2 (Convergence to a time-homogenous branching process) Assume that 
C{X_ n ) ==$■ fj, for some probability law fi on Ai ([0,1]). 

(a) IfO < 7* < 00, then 

c(x^ n , *_ n+ i,...) c(y£,yf ,...), (3.12) 

n— >oo 

where y^* is the time-homogenous branching process with log-Laplace operator Li^* in each 
step and initial law C(y$ ) = fi. 

(b) 7/7* = 0, then 

£ ((^(*)) t >o)„^c^(ra t >o)' ( 3 - 13 ) 

where denotes weak convergence of laws on path space, k n (t) := min{£; : < k < n, 
Yl?=k Ti — t}, an d y° is the super- Wright- Fisher diffusion with activity and growth parameter 
both identically 1 and initial law C(y^) = \i. 



The super- Wright-Fisher diffusion was studied in |FSfl3| . By definition, y° is the time-homo- 
geneous Markov process in Ai[0, 1] with continuous sample paths, whose Laplace functionals 
are given by 

E^[ e -( y tJ}] =e -(^Mtf) (/i G «M[0, 1], f £ B+[0, 1], t > 0). (3.14) 
Here IA® f = ut is the unique mild solution of the semilinear Cauchy equation 

I §- t u t (x) = \x(l - x)£zu t (x) + u t (x)(l - u t (x)) (t > 0, x £ [0, 1]), ^ ^ 

For a further study of the renormalization branching process X and its limiting processes y 7 
(7* > 0) we will use the technique of embedded particle systems, which we explain in the next 
section. 
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3.4 Weighted and Poissonized branching processes 

In this section, we explain how from a Poisson-cluster branching process it is possible to con- 
struct other branching processes by weighting and Poissonization. We first need to introduce 
spatial branching particle systems in some generality. 

Let E again be separable, locally compact, and metrizable. For v G N{E) and / G 
B[ i](E), we adopt the notation 

m m 

f°:=l and f u := JJ when v = Y, S *i (m > 1). (3.16) 

i=i i=i 

We call a continuous map x \— > Q(#, • ) from 12 into Mi(N(E)) a continuous offspring mech- 
anism. 

Fix continuous offspring mechanisms (1 < A; < n), and let (Xo, . . . ,X n ) be a Markov 
chain in M(E) such that, given that Xfc_i = Y^ILi the next step of the chain X& is a sum 
of independent random variables with laws Qk(xi, • ) (i = 1, • • • , ?7t). Then 

= (l _ C7i o...oU n f) u [y G M(E), f€B M (E)), (3.17) 
where £/& : 5r 0] i](-E) — > Sr 0j i](£?) is defined as 

t/ fc /(ar) := 1 - / Q fc (x,dt/)(1 - /)" (1 < A; < n, x G £, / G %i](£)). (3.18) 

We call Uk the generating operator of the transition law from X^-i to X&, and we call X = 
(Xq, . . . , X n ) the branching particle system on 2£ with generating operators Ui, . . . ,U n . It is 
often useful to write (|3.17|) in the suggestive form 

P v [Thm f (X n ) = 0] = P[Thin f/l0 ... oCW H = 0] (y G M(E), f G B [0 ,i](£))- (3.19) 

Here, if ^ is an jV(i?)-valued random variable and / G B^q^(E), then Thiny(i^) denotes an 
A^E 1 )- valued random variable such that conditioned on v, Thinj(z^) is obtained from v by 
independently throwing away particles from u, where a particle at x is kept with probability 
f(x). One has the elementary relations 

Tiu3Xf(Tiung(u)) = Thin /ff (i/) and Thi n/ (Pois(/i)) = Pois(//i), (3.20) 

where = denotes equality in distribution. 

We are now ready to describe weighted and Poissonized branching processes. Let X = 
(Xq, . . . ,X n ) be a Poisson-cluster branching process on E, with continuous weight functions 
qi, . . . ,q n , continuous cluster mechanisms Qi, ■ ■ ■ , Q n , and log-Laplace operators Ui, . . . ,U n 
given by ()3.2j) and satisfying (|3.3|) . Let Z% denote an A^(£')-valued random variable with 
law Q k (x, •)■ Let h G C+(E) be bounded, h ^ 0, and put E h := {x G E : h(x) > 0}. For 
/ G B + (E h ), define hf G B+(E) by hf(x) := h{x)f(x) if x G E h and hf(x) := otherwise. 

Proposition 3.3 (Weighting of Poisson-cluster branching processes) Assume that 
there exists a constant K < oo such that li^h < Kh for all k = 1, ... ,n. Then there exists a 
Poisson-cluster branching process X h = (Xq, . . . , X%) on E h with weight functions (q^, . . . , q%) 
given by q\ := qk/h, continuous cluster mechanisms Q\, . . . , Q 1 ^ given by 

Q h k (x,-):=£(hZ k x ) (x£E h ), (3.21) 
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and log-Laplace operators U^, . . . ,U% satisfying 

hU%f ■= MM) (f G B + {E h )). (3.22) 

The processes X and X h are related by 

C(Xft) = C(hX ) implies C(X%) = C(hX k ) (0 < k < n). (3.23) 

Proposition 3.4 (Poissonization of Poisson-cluster branching processes) Assume 
that U k h < h for all k = l,...,n. Then there exists a branching particle system X h = 
(Xq , . . . , X%) on E h with continuous offspring mechanisms Q\, . . . , Q\ given by 

Qh ^ - ):= m p[PoHhzkx) £ ' ] + (* ~ ffy)^ {x£ E% (3 - 24) 

and generating operators U± , . . . , U% satisfying 

hU%f :=U k (hf) (/ G B m (E h )). (3.25) 

The processes X and X h are related by 

C(X%) = £(Pois(hX )) implies C(X%) = C(Pois(hX k )) (0<k<n). (3.26) 

Here, the right-hand side of (|3.24|) is always a probability measure, despite that it may happen 
that q k {x)/h{x) > 1. The (straightforward) proofs of Propositions l3~3"l and l3~H can be found in 
Section [7.11 below. If (|3.23|) holds then we say that X h is obtained from X by weighting with 
density h. If ()3.26|) holds then we say that X h is obtained from X by Poissonization with 
density h. Proposition 13.41 savs that a Poisson-cluster branching process X contains, in a way, 
certain 'embedded' branching particle systems X h . Poissonization relations for superprocesses 
and embedded particle systems have enjoyed considerable attention, see |PSQ4j and references 
therein. 

A function h G B + {E) such that U k h < h is called lA k -superharmonic. If the reverse 
inequality holds we say that h is U k -subharmonic. \ilA k h = h then h is called li k -harmonic. 



3.5 Extinction versus unbounded growth for embedded particle systems 

In this section we explain how embedded particle systems can be used to prove Theorem 11.41 
Throughout this section ('jk)k>o are positive constants such that ^ n 7„ = oo and -y n — > 7* for 
some 7* G [0, 00), and X = (X_ n , . . . , Xq) is the renormalization branching process on [0, 1] 
defined in Section ET21 We write 

U^:=U ln _ 1 o...oU 10 . (3.27) 
In view of formula (|3.11j) . in order to prove Theorem 11.41 we need the following result. 
Proposition 3.5 (Limits of iterated log-Laplace operators) Uniformly on [0,1], 

(i) limW (n) p=l ( P eHi,i), 

(ii) limWWp = (p€Wo,o), (3.28) 

(iii) lim^)p = ^i 7 . (p€Wo,l), 

where pg,i 7 * : [0, 1] — > [0, 1] is a function depending on 7* but not on p G Ho,i- 
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In our proof of Proposition 13.51 we will use embedded particle systems X h = (X_ n , . . . ,X(j) 
obtained from X by Poissonization with certain h taken from the classes 7Yi 5 i, Wo,o> an d Wo,i- 

Lemma 3.6 (Embedded particle system with foi,i) The constant function h\^(x) := 1 
is -harmonic for each 7 > 0. The corresponding embedded particle system X 111 - 1 on [0, 1] 
satisfies 

P -M.[|xJm| g .] => 5^ (3.29) 

uniformly 2 for all x G [0, 1]. 

In Q3.29JI and similar formulas below, =>- denotes weak convergence of probability measures 
on [0, 00]. Thus, (|3.29|) says that for processes started with one particle on the position x at 
times — n, the number of particles at time zero converges to infinity as n — > 00. 

Lemma 3.7 (Embedded particle system with /io,o) The function ho t o(x) := x(l — x) 
(x G [0, 1]) is U-y-superharmonic for each 7 > 0. The corresponding embedded particle system 
X h °'° on (0, 1) is critical and satisfies 

p- n 'M|xJ°'°| G -1 5 Q (3.30) 
locally uniformly for all x G (0, 1). 

Here, a branching particle system X is called critical if each particle produces on average 
one offspring (in each time step and independent of its position). Formula (|3.3U|) says that 
the embedded particle system X h °'° gets extinct during the time interval {— n, ... ,0} with 
probability tending to one as n — * 00. We can summarize Lemmas 13 . 61 and 13 . 71 bv saying that 
the embedded particle system associated with grows unboundedly while the embedded 
particle system associated with /ir^o becomes extinct asn-> 00. 

We will also consider an embedded particle systems X 110 ' 1 for a certain /to,l taken from 
TCo,i- It turns out that this system either gets extinct or grows unboundedly, each with a 
positive probability. In order to determine these probabilities, we need to consider embedded 
particle systems for the time- homogeneous processes y 1 (7* G [0, 00)) from (|3.12|) and (|3.13() . 
If h G Wo,i is W 7 *-superharmonic for some 7* > 0, then Poissonizing the process y~< with 
h yields a branching particle system on (0, 1] which we denote by Y^* ,h = (Y 7 ' h , ' h , . . .). 
Likewise, if h G Wo,i is twice continuously differentiable and satisfies 

±x(l - x)j^h(x) - h(x)(l - h(x)) < 0, (3.31) 

then Poissonizing the super- Wright-Fisher diffusion with h yields a continuous-time branch- 
ing particle system on (0, 1], which we denote by Y 0,h = (Y® ,h )t>o- For example, for m > 4, 
the function h(x) := 1 — (1 — x) m satisfies (|3.31|) . 

Lemma 3.8 (Embedded particle system with /io,i) The function ho : i(x) := 1 — (1 — x) 7 
is Uj-superharmonic for each 7 > 0. The corresponding embedded particle system X 110 ' 1 on 
(0, 1] satisfies 

p -nA[| X Jo,i| G .] p r (x)8 OQ + (1- Pr (x))5 , (3.32) 

2 Since Ali [0,oo] is compact in the topology of weak convergence, there is a unique uniform structure 
compatible with the topology, and therefore it makes sense to talk about uniform convergence of 7Mi[0, 00]- 
valued functions (in this case, x t— > p~ n * Sx [|Xq 1,1 1 G ■ ] ). 
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locally uniformly for all x G (0, 1], where 

J P S X[Y V*M A ^ Q y k y 0] (0 < 7 * < oo)> 

Ay*(aO:=< (3.33) 

I p^[y t 0AM vt> o] (7* = o). 

We now explain how Lemmas l3.6H3.8l imply Proposition 13.51 In doing so, it will be more con- 
venient to work with weighted branching processes than with Poissonized branching processes. 
A little argument (which can be found in Lemma 17.121 below) shows that Lemmas l3.6H3.8l are 
equivalent to the next proposition. 

Proposition 3.9 (Extinction versus unbounded growth) Let hn, hoo, and /io,i be as 
in Lemmas \S. (MH. #1 For 7* G [0, 00) ; put p\ lri *{x) := 1, _Po,o,7*( x ) := (2; G [0,1]), and 

Po,i,7*(°) := and Pq,i,y( x ) '■= ho,i(x)p y *{x) (xG(0,l]), (3.34) 
with p 7 » as in kS.Sty . Then, for (I, r) = (1, 1), (0, 0), and (0, 1), 

p-^[(^o,V)e-] e~ p trA x h + {1-6-^,1*^)5^ (3.35) 

L J n— »oo 

uniformly for all x G [0, 1] . 

Formula (|3.35|) says that the weighted branching process X hl ' r exhibits a form of 'extinction 
versus unbounded growth'. More precisely, for large n the total mass of Hi t Xq is close to or 
00 with high probability. 

Proof of Proposition iniU By (l3~4l . 

U^p(x) = -logE- n ' d *[e-( X 0>P}] (p G 5+ [0,1], x G [0,1]). (3.36) 
We first prove formula 1)3. 28 j) (ii). For (Z, r) = (0,0), formula (|3.35|) says that 

P~ n ' 5x [(X ,h 00 ) G •] =>- 5 (3.37) 

n—*oo 

uniformly for all x G [0,1]. If p G TCofl, then we can find r > such that p < rhofl. Therefore, 
(|3.37j) implies that for any p G Hq,o, 



j— n,<5aD 

By (|3~36|) it follows that 



p- n < d *[{X ,p) G •] So. (3.38) 

n— >oo 



U { - n) p(x) = -\ogE- n ^\e-^P)} — ► 0, (3.39) 

where the limits in (|3.38j) and (|3.39|) are uniform in x G [0, 1]. This proves formula (|3.28j) (ii). 
To prove formula (|3.28|) (in), note that for any p G Ho,l we can choose < r_ < r + such that 
f-ho,i < P + ^0,0 < r +^o,i- Therefore, (|3.35f) implies that 

p- n ^[(X ,p) + (X ,h 0fi ) €■] => e -Po,i l7 *(^) 5o + (i_ e-Po.i.r^))^. (3.40) 
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Using moreover (j3.37|) . we see that 

p- n > s *[(X ,p) G •] => e~ p *°> l >r( x h + (1 - e^Lr^Voo- (3.41) 

n— >oo 

By (l3~3fi|l . it follows that 

U^p(x) = -logE~ n ^\e-^,ph — , p * f x \ (3.42) 

where all limits are uniform in x G [0, 1]. This proves (|3.28f) (iii). The proof of (|3.28j) (i) is 
similar but easier. I 

4 Discussion, open problems 

4.1 Discussion 

Consider a ([0, 1] 2 ) z2 -valued process x = (x^)^ gZ 2 = (x^,x|)^ gZ 2, solving a system of SDE's 
of the form 



dx|(i) = Yl (xi(t)-xi(t))dt+^2axi(t)(l-xi(t))d^(t), 
v- l»?-£|=i 



dx|(i)= £ (x 2 (i)-x|(t))di+^2p(xi(f))x|(i)(l-x|(t))di?|(t), 
r?:[77-e|=l 



(4.1) 



where a > is a constant, p is a nonnegative function on [0, 1] satisfying p(0) = and p(l) > 0, 
and (-B|)^~22 is a collection of independent Brownian motions. We call x a system of linearly 
interacting catalytic Wright-Fisher diffusions with catalyzation function p. It is expected that 
x clusters, i.e., x(i) converges in distribution as t — > oo to a limit (x^(oo))^ g ^2 such that 
x^(cxd) = xo(oo) for all £ G Z 2 and xo(oo) takes values in the effective boundary associated 
with the diffusion matrix w a ' p (see (|2.3|1 ). Heuristic arguments, based on renormalization, yield 
a formula for the clustering distribution £(xo(oo)) in terms of the diffusion matrix w* which 
is the unique solution of the asymptotic fixed point equation (J2.16|) (ii) in the renormalization 
class W^t! see Conjecture IA.3I in Appendix IA. 21 below. 

The present paper is inspired by the work of Greven, Klenke and Wakolbinger GK WOl] . 
They study a model that is closely related to (|4.1[) . but where x 1 is replaced by a voter model. 
They show that their model clusters and determine its clustering distribution £(xq(oo)), which 
turns out to coincide with the mentioned prediction for (j4.1|) based on renormalization theory. 
In fact, they believe their results to hold for the model in (|4,1|) too, but they could not prove 
this due to certain technical difficulties that a [0, l]-valued catalyst would create, compared 
to the simpler {0, l}-valued voter model. 

The work in jCKWOll not only provides the main motivation for the present paper, but 
also inspired some of our techniques for proving Theorem 11.41 This concerns in particular the 
proof of Proposition l3.11 which makes the connection between renormalization transformations 
and a branching process. We hope that conversely, our techniques may shed some light on the 
problems left open by |(tKW01| . in particular, the question whether their results stay true if 
the voter model catalyst is replaced by a Wright-Fisher catalyst. It seems plausible that their 
results may not hold for the model in (j4.1j) if the catalyzing function p grows too fast at 0. On 
the other hand, our proofs suggest that p with a finite slope at should be OK. (In particular, 
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while deriving formula (|3.40|) . we use that p can be bounded from above by r+/io,l f° r some 
r + > 0, which requires that p has a finite slope at 0.) 

Our results are also interesting in the wider program of studying renormalization classes 
in the sense of Definition 11.11 We conjecture that the class W^t, unlike all renormalization 
classes studied previously, contains no fixed shapes (see the discussion following Lemma l2.8[) . 
In fact, we expect this to be the usual situation. In this sense, the renormalization classes 
studied so far were all of a special type. 

4.2 Open problems 

The general program of studying renormalization classes in the sense of Definition contains 
a wealth of open problems. In our proofs, we make heavy use of the single-way nature of the 
catalyzation in (|1.7|) . in particular, the fact that y 1 is an autonomous process which allows 
one to condition on y 1 and consider y 2 as a process in a random environment created by y . 
As soon as one leaves the single-way catalytic regime one runs into several difficulties, both 
technically (it is hard to prove that a given class of matrices is a renormalization class in the 
sense of Definition II. 1J) and conceptually (it is not clear when solutions to the asymptotic 
fixed shape equation (|2.16|) (ii) are unique). Therefore, it seems at present hard to verify the 
complete picture for renormalization classes on the unit square that arises from the numerical 
simulations described in Section 12.21 and Figures [21 and unless one or more essential new 
ideas are added. 

In this context, the study of the nonlinear partial differential equation (j2.18|) and its fixed 
points seems to be a challenging problem. This may be a hard problem from an analytic point 
of view, since the equation is degenerate and not in divergence form. For the renormalization 
class Wcatj the quasilinear equation (|2.18|) reduces to the semilinear equation (|3.15|) . which 
is analytically easier to treat and moreover has a probabilistic interpretation in terms of a 
superprocess. For a study of the semilinear equation (|3.15|) we refer to |FS03| . We do not 
know whether solutions to equation (|2.18j) can in general be represented in terms of a stochastic 
process of some sort. 

Even for the renormalization class W ca t) several interesting problems are left open. One of 
the most urgent ones is to prove that the functions p$ 1 7 * are not constant in 7*, and therefore, 
by Lemma 12.81 (c). contains no fixed shapes. Moreover, we have not investigated the 

iterated renormalization transformations in the regime 7* = 00. Also, we believe that the 
convergence in (j3,28|) (ii) does not hold if the condition that p is Lipschitz is dropped, in 
particular, if p has an infinite slope at or an infinite negative slope at 1. For p £ TCo t o, it seems 
plausible that a properly rescaled version of the iterates U^p converges to a universal limit, 
but we have not investigated this either. Finally, we have not investigated the convergence of 
the iterated kernels K w >( n > from (|2.4j) (in particular, we have not verified Conjecture IA.2|) for 
the renormalization class W ca t- 

Our methods, combined with those in BCGdH95 , can probably be extended to study the 
action of iterated renormalization transformations on diffusion matrices of the following more 
general form (compared to (|l-4(l ): 

-W-^?' (*=€[<U1 2 ), (42) 

where g : [0, 1] — ► R is Lipschitz, g(0) = g(l) = 0, g > on (0,1), and p £ TL as before. 
This would, however, require a lot of extra technical work and probably not generate much 
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new insight. The numerical simulations mentioned in Section \'2.2\ suggest that many diffusion 
matrices of an even more general form than 1)4. 2|) also converge under renormalization to the 
limit points w* from Theorem 11.41 but we don't know how to prove this. 

Part II 

Outline of Part II In Section we verify that W C at is a renormalization class, we prove 
Proposition 13.11 which connects the renormalization transformations F c to the log-Laplace 
operators Z/y, and we collect a number of technical properties of the operators U-y that will be 
needed later on. In Section we prove Theorem 13.21 about the convergence of the renormal- 
ization branching process to a time-homogeneous limit. In Section we prove the statements 
from Section 13.51 about extinction versus unbounded growth of embedded particle systems, 
with the exception of Lemma 13.71 which is proved in Section |SJ In Section EE finally, we 
combine the results derived by that point to prove our main theorem. 

5 The renormalization class W ca t 

In this section we prove Theorem II .41 (a) and Proposition l3.il as well as Lemmas 12. 1H2.8I from 
Section[21 The section is organized according to the techniques used. Section 15 . II collects some 
facts that hold for general renormalization classes on compact sets. In Section [5.21 we use the 
SDE (|1.7j) to couple catalytic Wright-Fisher diffusions. In Section 15.31 we apply the moment 
duality for the Wright-Fisher diffusion to the catalyst and to the reactant conditioned on the 
catalyst. In Section ^. 4l we prove that monotone concave catalyzing functions form a preserved 
class under renormalization. 

5.1 Renormalization classes on compact sets 

In this section, we prove the lemmas stated in Sectional Recall that D C M rf is open, bounded, 
and convex, and that W is a prerenormalization class on D, equipped with the topology of 
uniform convergence. 

Proof of Lemma l2.1I To see that (x, c, w) i— > u x w is continuous, let (x n , c n , w n ) be a sequence 
converging in D x (0, oo) x W to a limit (x,c,w). By the compactness of D, the sequence 
(vxZ' W ™)n>o is tight, and each limit point v* satisfies 

(v*,AZ™f) = (f<=C®(D)). (5.1) 

Therefore, by |KK86I Theorem 4.9.17], v* is an invariant law for the martingale problem 
associated with A x w ■ Since we are assuming uniqueness of the invariant law, v* = v x and 
therefore v x n n ,Wn =^ v< x" • The continuity of F c w(x) is a simple consequence of the continuity 

r C,W _ 

ot u x ' . I 

Proof of Lemma 12.21 Formula H2.1|) (i) follows from the fact that rescaling the time in 
solutions (y*)t>o to the martingale problem for A x w by a factor A has no influence on the 
invariant law. Formula (|2.1jl (ii) is a direct consequence of formula (|2.1|) (i). I 

Proof of Lemma 12. 3I This follows by inserting the functions f(x) = X{ and fix) = XiXj into 
the equilibrium equation ()5.1|) . I 

Proof of Lemma 12.41 If x £ d w D, then yt := x (t > 0) is a stationary solution to the 
martingale problem for A x ' w , and therefore v x w = S x and F c w(x) = w(x) = 0. On the other 
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hand, if x ^ d w D, then yt := x (t > 0) is not a stationary solution to the martingale problem 
for Ax' w and therefore fjyVx™ (dy)\y — x\ 2 > 0. Let tr(u>(y)) := Yli w a(y) denote the trace 
of w(y). By Q (ii), ifr^)^) = |/^i/^(dy)tr(«;(y)) = /^^(dy)^ - x| 2 > and 
therefore F c w(x) ^ 0. I 

From now on assume that W is a renormalization class. Note that 

K w,(n) =J/ c n - 1 ,F^-Vw... J/ c ,w (n>l), (5.2) 

where we denote the composition of two probability kernels K, L on D by 

(KL) x (dz) := (_K x (dy)L y (dz). (5.3) 
Jd 

Proof of Lemma 12.51 This is a direct consequence of Lemmas 12.11 and 12.31 In particular, the 

relations (|2,6|) follow by iterating the relations (|2.2j) . I 

Proof of Lemma 12.61 Recall that ti{w{y)) denotes the trace of w{y). Formulas (|2.5j) and 
((23)) (ii) show that 

iC' (n) (dy) |y ~x\ 2 = s n f_K^ n \dy)tT{w{y)). (5.4) 



D 



Since D is compact, the left-hand side of this equation is bounded uniformly in x G D and 
n > 1, and therefore, since we are assuming s n — > oo, 

lim sup / ^'^(dy)tr(w(y)) = 0. (5.5) 

Since u> is symmetric and nonnegative definite, tr(w(y)) is nonnegative, and zero if and only if 
y G d w D. If / G C(D) satisfies / = on <9„,-D, then, for every e > 0, the sets C m := {x £ D : 
\f(x)\ > e + mti(w(x))} are compact with C m j as m | oo, so there exists an m (depending 
on e) such that |/| < e + mtr(w). Therefore, 



lim sup sup 



K^ n \dy)f{y) <limsupsup _K^ n \dy)\f(y)\ 



l->oo xeD Jfl ra->oo x€ dJD 

< e + mlimsup sup j_K x ^ n \dy)ti(w(y)) = e. 



(5.6) 



Since e > is arbitrary, (|2.7|) follows. I 

Proof of Lemma Hm By (l2~TUl) . (t2~T2l . and (|2~T3f) . w** = lim^oo^*)™^ for each w G W. 
By Lemma l2~T1 fb). F 7 * : W — > W is continuous, so w** is the unique fixed point of F 7 *. This 
proves part (a). 

Now let 7^ w G W and assume that W = {Xw : A > 0} is a fixed shape. Then 
W 3 s n F^ w — > w** whenever s n — > oo and s n+ i/s n -^1 + 7* for some < 7* < 00, which 

shows that VV = {Xw*, : X > 0}. Thus, W can contain at most one fixed shape, and if it does, 
then the w** for different values of 7* must be constant multiples of each other. This proves 
part (c) and the uniqueness statement in part (b). 

To complete the proof of part (b), note that if w* = w** does not depend on 7*, then 
w* G W solves ()2.16j) (i) for all < 7* < 00, hence F c w* = (1 + \)~ l w* for all c > 0, and 



therefore, by scaling fLemmal2~2l. F c (Xw*) = XF c/x (w*) = A(l + ^)~ 1 w* = (j + lw 



w 
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5.2 Coupling of catalytic Wright-Fisher diffusions 

In this section we verify condition (i) of Definition 11.11 for the class W ca t , and we prepare for 
the verification of conditions (ii)-(iv) in Section 15.31 In fact, we will show that the larger 
class W ca t := {w a ' p : a > 0, p G C+[0, 1]} is also a renormalization class, and the equivalents 
of Theorem 11.41 (a) and Proposition 13.11 remain true for this larger class. (We do not know, 
however, if the convergence statements in Theorem 11.41 (b) also hold in this larger class; see 
the discussion in Section POD 

For each c > 0, w G W cat and x G [0, l] 2 , the operator A^ w is a densely defined linear 
operator on C([0, l] 2 ) that maps the identity function into zero and, as one easily verifies, 
satisfies the positive maximum principle. Since [0, l] 2 is compact, the existence of a solution 
to the martingale problem for A% w , for each [0, l] 2 -valued initial condition, now follows from 
general theory (see |RW87j . Theorem 5.23.5, or |EK86| Theorem 4.5.4 and Remark 4.5.5]). 

We are therefore left with the task of verifying uniqueness of solutions to the martingale 
problem for Ax' w . By EK86, Problem 4.19, Corollary 5.3.4, and Theorem 5.3.6], it suffices to 
show that solutions to (|1,7[) are pathwise unique. 

Lemma 5.1 (Monotone coupling of Wright-Fisher diffusions) Assume that < x < 
x < 1, c > and that (Pt)t>o is a progressively measurable, nonnegative process such that 
su Pt>o wen Pt{w) < oo. Let y,y be [0, l]-valued solutions to the SDE's 



dy t = c (x - y t )dt + y/2P t y t {\ - YtjdB t , ( 
dy t = c(z - y t )dt + y/2P t y t (l - y t )dB t , 

where in both equations B is the same Brownian motion. If yo < yo a -S-, then 

yt<yt Vt>0 a.s. (5.8) 

Proof This is an easy adaptation of a technique due to Yamada and Watanabe |YW71j . Since 
Jo+ = °°' ^ ^ s P oss ibl e to choose p n G C[0, oo) such that J °° p n (x)dx = 1 and 

< p n (x) < —1(0,1] (x) (x > 0). (5.9) 
nx 

Define <j> n G C^(R) by 

<j) n {x) := / dy dz Pn {z). (5.10) 



JO J0 

One easily verifies that (j> n (x), x(p' n (x), and xcpn(x) are nonnegative and converge, as n — > oo, 
to x V 0, x V 0, and 0, respectively. By Ito's formula: 

E[Myt - yt)]=E[Myo - yo)] t (i) 

+c (x-x) I E\4>' n {y s - y s )}ds -c E[(y s - y s )^(y s ~ ys)}ds (ii) 
t Jo Jo 

+ J E[p s y Ys (i-y s )- Vy s (i-y s )) 2 <(y 8 -y s )]ds. (m) 

(5.11) 

Here the terms in (ii) are nonpositive, and hence, letting n — * oo and using the elementary 
estimate 

Wy(i-y)-Vy(i-y)\ < (y,y g [0,1]), (5.12) 
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the properties of <f) n , and the fact that the process P is uniformly bounded, we find that 

E[0 V (y ( - y t )] < E[0 V (y - yo)] = 0, (5.13) 

by our assumption that yo < yo- This shows that yt < yt a.s. for each fixed t > 0, and by 
the continuity of sample paths the statement holds for all t > almost surely. I 

Corollary 5.2 (Pathwise uniqueness) For all c > 0, a > 0, p G C + [0, 1] and x G [0,1], 
solutions to the SDE fl 1. 7| ) are pathwise unique. 

Proof Let (y^y 2 ) and (y 1 ,^ 2 ) be solutions to Q1-7JI relative to the same pair (B 1 ,^ 2 ) of 
Brownian motions, with (y^y 2 ,) = (yo>yo)- Appfyi n g Lemma I5.H with inequality in both 
directions, we see that y 1 = y 1 a.s. Applying Lemma l5.1l two more times, this time using that 
y 1 = y 1 a.s., we see that also y 2 = y 2 a.s. I 

Corollary 5.3 (Exponential coupling) Assume that x G [0, 1], c > 0, and a > 0. Let y,y 

be solutions to the SDE 



dy t = c{x- y t )dt + v / 2ay t (1 - y t )dB t , (5.14) 

relative to the same Brownian motion B. Then 

E[\yt-yt\] =e- ct E[\y -y \}. (5.15) 

Proof If yo = y and yo = y are deterministic and y < y, then by Lemma 15.11 and a simple 
moment calculation 

E[\yt - yt\] = E[y t - y t ] = e~ ct \y - y\. (5.16) 

The same argument applies when y > y. The general case where yo and yo are random follows 
by conditioning on (yo,yo)- ■ 

Corollary 5.4 (Ergodicity) The Markov process defined by the SDE I.V. 6|) has a unique 
invariant law Tj and is ergodic, i.e, solutions to h'J. b]) started in an arbitrary initial law £(yo) 
satisfy £(y t ) => Tj. 

Proof Since our process is a Feller diffusion on a compactum, the existence of an invariant 
law follows from a simple time averaging argument. Now start one solution y of (|3.6|) in 
this invariant law and let y be any other solution, relative to the same Brownian motion. 
Corollary 15.. 31 then gives ergodicity and, in particular, uniqueness of the invariant law. I 

Remark 5.5 (Density of invariant law) It is well-known (see, for example |Ewe04| for- 
mula (5.70)]) that TZ is a f3(a%, «2)-distribution, where a\ := x/7 and 02 := (1 — 2O/7, i-e., 
Vl = 5 X (x G {0, 1}) and 

TZ{dy) = rfclnal) yai ~ 1{1 " yr ~ ldy {x G (0 ' (5 ' 17) 

o 
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We conclude this section with a lemma that prepares for the verification of condition (iv) in 
Definition 11.11 for the class W ca t- 

Lemma 5.6 (Monotone coupling of stationary Wright-Fisher diffusions) Assume 
that c > 0, a > and < x < x < 1. Then the pair of equations 

dy t = c(x - y t )dt + y / 2ay t (l - y t )dB t , _ 
dy t = c(x- y t )dt + y/2ay t (l - y t )dB t 

has a unique stationary solution (yt,yt)teR- This stationary solution satisfies 

yt<Yt VteR a.s. (5.19) 

Proof Let (yt,yt)t>o be a solution of (|5.18l) and let (y' t ,y' t )t>o be another one, relative to the 
same Brownian motion B. Then, by Lemma l5.Hl E[\yt — y' t \] — > and also E[\yt — y' t \] — > as 
t —* oo. Hence we may argue as in the proof of Corollary 15 . 41 that ()5.18(l has a unique invariant 
law and is ergodic. Now start a solution of (|5.18j) in an initial condition such that yo < yo- 
By ergodicity, the law of this solution converges as t — ► oo to the invariant law of (|5.18|) and 
using Lemma 15.11 we see that this invariant law is concentrated on {(y,y) 6 [0, l] 2 : y < y}. 
Now consider, on the whole real time axis, the stationary solution to (|5.18|) with this invariant 
law. Applying Lemma 15. II once more, we see that (EHflj) holds. I 



5.3 Duality for catalytic Wright-Fisher diffusions 

In this section we prove Theorem 11.41 (a) and Proposition 13.11 Moreover, we will show that 
their statements remain true if the renormalization class W C at is replaced by the larger class 
VVcat := {w a,p : a > 0, p E C+[0, 1]}. We begin by recalling the usual moment duality for 
Wright-Fisher diffusions. 

For 7 > and x £ [0, 1], let y be a solution to the SDE 

dy(t) = i (x - y(t))dt + v/2y(t)(l " y(t))dB(t), (5.20) 

i.e., y is a Wright-Fisher diffusion with a linear drift towards x. It is well-known that y has 
a moment dual. To be precise, let (<fi, ip) be a Markov process in N 2 = {0, 1, . . .} 2 that jumps 
as: 

{4>t,i>t) (4>t ~ 1, tpt) with rate <j) t (<Pt - 1) (f . „ n 

(0 t ,Vt)^(^ -1,^ + 1) with rate 

Then one has the following duality relation (see for example Lemma 2.3 in |Shi80| or Propo- 
sition 1.5 in IGKWOlj ) 

EV [y?x m ] = S (n ' m) [y^x^] (y G [0, 1], (n, m) G N 2 ), (5.22) 

where 0° := 1. The duality in (|5.22|) has the following heuristic explanation. Consider a 
population containing a fixed, large number of organisms, that come in two genetic types, 
say I and II. Each pair of organisms in the population is resampled with rate 2. This means 
that one organism of the pair (chosen at random) dies, while the other organism produces one 
child of its own genetic type. Moreover, each organism is replaced with rate i by an organism 
chosen from an infinite reservoir where the frequency of type I has the fixed value x. In the 
limit that the number of organisms in the population is large, the relative frequency yt of type 
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I organisms follows the SDE (|5.2U|) . Now E[y™] is the probability that n organisms sampled 
from the population at time t are all of type I. In order to find this probability, we follow 
the ancestors of these organisms back in time. Viewed backwards in time, these ancestors 
live for a while in the population, until, with rate ^, they jump to the infinite reservoir. 
Moreover, due to resampling, each pair of ancestors coalesces with rate 2 to one common 
ancestor. Denoting the number of ancestors that lived at time t — s in the population and in 
the reservoir by (f) s and ip s , respectively, we see that the probability that all ancestors are of 
type I is Ey[y™\ = E^[y^x^}. This gives a heuristic explanation of ()5.22|) . 

Since eventually all ancestors of the process (</>, ijj) end up in the reservoir, we have 
{(j)t,ipt) (0,^00) as t ~~ * 00 a - s - f° r some N-valued random variable Voo- Taking the limit 
t — > 00 in ([5.22ft . we see that the moments of the invariant law Tj from Corollary 15 .41 are given 
by: 

' TZ(dy)y n = £ ( "' 0) [x*°°] (n > 0). (5.23) 



It is not hard to obtain an inductive formula for the moments of T%, which can then be solved 
to yield the formula 

I Tl{dy)y n = J] (n > 1). (5.24) 



In particular, it follows that 



/ 



TZ(dy)y{l-y) = -—x{l-x). (5.25) 
1 + 7 



This is the important fixed shape property of the Wright-Fisher diffusion (see formula 112.171) ). 

We now consider catalytic Wright-Fisher diffusions (y x ,y 2 ) as in (|1.7|) with p G C+[0, 1] 
and apply duality to the catalyst y 2 conditioned on the reactant y 1 . Let (yl ,yt)teR be a 
stationary solution to the SDE (|1,7|) with c = I/7. Let (4>,ip) be a N 2 -valued process, defined 
on the same probability space as (y^y 2 ), such that conditioned on the past path (yl t )t<o, 
the process (<p, ip) is a (time-inhomogeneous) Markov process that jumps as: 



$j,i>j) i<th ~ hih) witn rate PiY-MMt - 1), 

(4>t,ipt) -> {4>t - 1, $t + 1) with rate ±<j) t . 



(5.26) 



Then, in analogy with (|5.22[) . 

^[(yo) n ^l(y- f )*<o] = E^[{yl t )^xf\{yl t ) t < Q ] ((n,m) G N 2 , t > 0). (5.27) 

We may interpret (|5.26f) by saying that pairs of ancestors in a finite population coalesce with 
time-dependent rate 2p(y! t ) and ancestors jump to an infinite reservoir with constant rate 
-. Again, eventualy all ancestors end up in the reservoir, and therefore (^>ti^t) ~^ (0, ?/>oo) as 
t — > 00 a.s. for some N-valued random variable ipoo. Taking the limit t — » 00 in (|5.27|) we find 
that 

£[(yo) n ^l(y-t)*<o] = S (n ' m) [4°°l(y-*)t<o] ((n,m) G N 2 , t > 0). (5.28) 



Lemma 5.7 (Uniqueness of invariant law) For each c > 0, w G W ca t, and x G [0, l] z 

there exists a unique invariant law Vx' w for the martingale problem for Ax' w . 
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Proof Our process being a Feller diffusion on a compactum, the existence of an invariant 
law follows from time averaging. We need to show uniqueness. If (y^y 2 ) = yl,y 2 )t&R is a 
stationary solution, then y 1 is an autonomous process, and C(yq) = T l J c the unique invariant 
law from Corollary 15.41 Therefore, C((yj )t£m) is determined uniquely by the requirement 
that (y^y 2 ) be stationary. By (|5.28j) . the conditional distribution of y 2 , given (yl)t<o is 
determined uniquely, and therefore the joint distribution of y 2 , and (y^ )«o is determined 
uniquely. In particular, £(yg,yg) = v% w is determined uniquely. I 



Remark 5.8 (Reversibility) It seems that the invariant law Vx W from Lemma 15.71 is re- 
versible. In many cases (densities of) reversible invariant measures can be obtained in closed 
form by solving the equations of detailed balance. This is the case, for example, for the one- 
dimensional Wright-Fisher diffusion. We have not attempted this for the catalytic Wright- 
Fisher diffusion. O 

The next proposition implies Proposition 13.11 and prepares for the proof of Theorem 11.41 (a). 

Proposition 5.9 (Extended renormalization class) The set Wcat is a renormalization 
class on [0, l] 2 , and 

F 1 w l ^ p = w 1 ^ p (p€C+ [0,1], 7>0). (5.29) 

Proof To see that W ca t is a renormalization class we need to check conditions (i)-(iv) from 
Definition 11.11 By Lemma 15.21 the martingale problem for A x ' w is well-posed for all c > 0, 
w G Wcat and x G [0, l] 2 . By Lemma f5.71 the corresponding Feller process on [0,1] 2 has 
a unique invariant law Vx ■ This shows that conditions (i) and (ii) from Definition 11.11 are 
satisfied. Note that by the compactness of [0, l] 2 , any continuous function on [0, l] 2 is bounded, 
so condition (iii) is automatically satisfied. Hence W is a prerenormalization class. As a 
consequence, for any p 6 C+[0, 1], F~ 1 w 1,p is well-defined by (jl.2|) and (|2,9[) . We will now first 
prove (|5.29|) and then show that W ca t is a renormalization class. 

Fix 7 > 0, p G C+[0, 1], and x G [0, l] 2 . Let (yj,y 2 )j g R be a stationary solution to the 
SDE Cd) with a = 1 and c = I/7. Then 

F,wlf(x) = (l+7)E[wlf(ylyl)] (ij = 1,2). (5.30) 

Since = if % ^ j, it is clear that Fjw]j P (x) = if i ^ j. Since £(yj) = it follows 
from (l5~231) that F^w{f (x) =x 1 (l-x 1 ). We are left with the task of showing that 



Here, by J22I) (ii), 



F^w^fix) = U 1 p{x\)x2{i — X2). (5.31) 
F 1 w^{x) = {l + 1 )E[p{yl)yl{l-yl)] 



;i + l)E[(y 2 - x 2 ] 



2l 



(5.32) 



By (|5.28D . using the fact that -E[yg] = £2 (which follows from (|5.27|) or more elementary from 
(TIoT) (i)), we find that 

E[(y 2 o -*2?\ = E[{yl?\ - {X2? = E<W\xf°) - (x 2 ) 2 = = l]x 2 (l-x 2 ) (t > 0). 

(5.33) 
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Note that p( 2 '°) [ip^ = 1] is the probability that the two ancestors coalesce before one of them 
leaves the population. The probability of noncoalescence is given by 



pWtfn = 2}= £[e" /<T 7 My-t)dty 



(5.34) 



where r 7 is an exponentially distributed random variable with mean 7. Combining this with 
(tQ2l and (l5~331l we find that 



Fy 2 f(x) = a + 1)E[1 - e & " Piy -t/2 )dt ]x 2 (l - x 2 ) 



--q 7 E[l - e~( Zx 
-Ujp(xi)x2(l - x 2 ), 



x 2 (l - x 2 ) 



(5.35) 



where we have used the definition of IA^. 

We still have to show that W cat satisfies condition (iv) from Definition ll.il For any a > 
and p e C+[Q, 1], by scaling (Xemma l2~2j) and (|5.29[) . 



™ « 1 £■ , a, 1 — 1 £ 

F c w a ' p = aFsLtv ' a = a(l + -)" i Fcu; 1 ' q 



(5.36) 



By Lemma l2.1| this diffusion matrix is continuous, which implies that is continuous.! 

Our proof of Propostion 15.91 has a corollary. 

Corollary 5.10 (Continuity in parameters) The map (a;, 7) 1— > Q 7 (av) /ram [0,1] x 
(0,oo) to .Mi(.M [0, 1]) and the map (x , 7 , p) h-> U^p(x) from [0, 1] x (0, 00) x C+[0, 1] to M are 
continuous. 

Proof By Lemma l2.1| the diffusion matrix in l)5.36jl is continuous in x, 7, and p, which implies 

the continuity of U 7 p(x). It follows that the map is continuous 

for all / G C_|_[0, 1], so by |Kal761 Theorem 4.2], (£,7) 1— > Q 7 (x, •) is continuous. I 

Proof of Theorem 11.41 (a) We need to show that W ca t is a renormalization class and that 



F c maps the subclasses W cat into themselves. It has already been explained in Section El that 
the latter fact is a consequence of Lemma 12.41 Since in Proposition 15.91 it has been shown 
that W ca t is a renormalization class, we are left with the task to show that F c maps W C at into 
itself. By (|5.29j) and scaling, it suffices to show that IA 7 maps 7i into itself. 

Fix < x < x < 1. By Lemma 15. 61 we can couple the processes y2 and y~ from (|3.6|) such 
that 

yZ(t)<yl(t) Vt<0 a.s. (5.37) 
Since the function z \—* 1 — e~ z on [0, 00) is Lipschitz continuous with Lipschitz constant 1, 
\Ujp(x) — U 1 p{x) J 

(I + 1)E [i _ e -/ Q T >(ylH/2))dt ] _ (i + 1)E[1 _ e - J? ptyl(-t/2))dt] 



\p(yl(-t/2)) - p(y2(-t/2))\dt 
7 |y?(-t/2)-y2(-t/2)|dt 



(5.38) 



= (I + l)Isy(x -x) = L(l + 7 )|z - x|, 

where L is the Lipschitz constant of p and we have used the same exponentially distributed 
r 7 for yl and y~. I 
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5.4 Monotone and concave catalyzing functions 

In this section we prove that the log-Laplace operators lA^ from (|3.9|) map monotone functions 
into monotone functions, and monotone concave functions into monotone concave functions. 
We do not know if in general maps concave functions into concave functions. 

Proposition 5.11 (Preservation of monotonicity and concavity) Let 7 > 0. Then: 

(a) If f G C+[0, 1] is nondecreasing, thenU^f is nondecreasing. 

(b) // / G C+[0, 1] is nondecreasing and concave, then lA^f is nondecreasing and concave. 

Proof Our proof of Proposition 15 . 1 ll is in part based on ideas from JBCGdH9fJ Appendix A]. 
The proof is quite long and will depend on several lemmas. We remark that part (a) can be 
proved in a more elementary way using Lemma 15.61 

We recall some facts from Hille-Yosida theory. A linear operator A on a Banach space V 
is closable and its closure A generates a strongly continuous contraction semigroup (St)t>o if 
and only if 

(i) T>(A) is dense, 

(ii) A is dissipative, (5.39) 

(iii) 1Z(1 — aA) is dense for some, and hence for all a > 0. 

Here, for any linear operator B on V, T>(B) and Tt{B) denote the domain and range of B, 
respectively. For each a > 0, the operator (1 — aA) : T>(A) — > V is a bijection and its inverse 
(1 — aA) -1 : V — ► T>(A) is a bounded linear operator, given by 

POD 

(1 - aA)- l u = \ S t u a~ l e- t/a dt (u G V, a> 0). (5.40) 
Jo 

If E is a compact metrizable space and C{E) is the Banach space of continuous real functions 
on E, equipped with the supremumnorm, then a linear operator A on C(E) is closable and its 
closure A generates a Feller semigroup if and only if (see [EK86, Theorem 4.2.2 and remarks 
on page 166]) 

(i) 1 G V(A) and ll = 0, 

(ii) V(A) is dense, , g ^ 

(iii) A satisfies the positive maximum principle, 

(iv) TZ(1 — aA) is dense for some, and hence for all a > 0. 

If A generates a Feller semigroup and g G C(E), then the operator A + g (with domain 
V(A + g) := T>(A)) generates a strongly continuous semigroup (Sf)t>o on C(E). If g < 
then (Sf)t>o is contractive. If (£t)t>o is the Feller process with generator A, then one has the 
Feynman-Kac representation 

Sfu(x) = E x lu(Z{t))e-fo 9 ^( s )) ds ] (t>0,xeE,g,ue C{E)). (5.42) 

Let C( n )([0,l] 2 ) denote the space of continuous real functions on [0, l] 2 whose partial deriva- 
tives up to n-th order exist and are continuous on [0, l] 2 (including the boundary), and 
put C(°°)([0, l] 2 ) := f|„C (n) ([0, 1] 2 )- Define a linear operator B on C([0,1] 2 ) with domain 
V(B) :=C(°°)([0,1] 2 ) by 

Bu(x, y) := y(l - y)^u(x, y) + ±{x - y)-§^u(x, y). (5.43) 

Below, we will prove: 
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Lemma 5.12 (Feller semigroup) The closure in C([0, l] 2 ) of the operator B generates a 
Feller semigroup on C([0, l] 2 ). 



Write 



C+ :={u G C([0,1] 2 ) : u > 0}, 
Cu 
C2^ 



{nECW([0,l] 2 ):|n,|n>0}, (5.44) 
{.EC( 2 )([0,l] 2 ):^n,^,^n>0}. 



Let 5 denote the closure of a set S C C([0, 1] ). We need the following lemma. 

Lemma 5.13 (Preserved classes) Let g G C([0, l] 2 ) and let (Sf)t>o be the strongly contin- 
uous semigroup with generator B + g. Then, for each t > 0: 

(a) If g G C\ + , then Sf maps C + n C\ + into itself. 

(b) // g G Ci + n C2+, then Sf maps C + n Ci+ fl C2+ into itself. 

To see why Lemma 15. 131 implies Proposition 15. 1 ll let (x(t),y(t))t>o denote the Feller process 
in [0, l] 2 generated by B. It is easy to see that x(i) = x(0) a.s. for all t > 0. For fixed x(0) = x, 
the process (y(t))t>o is the diffusion given by the SDE (|5.2U|) . Therefore, by Feynman-Kac, 
for each g G C([0,1] 2 ), 

^[ e jJ^,y(s))ds] = 5fl(x,y), (5.45) 
where 1 denotes the constant function 1 G C([0, l] 2 ). By (|3.9|) . 

^ 7 /(x) = (i + l)(l-|r2(dy)^[e-/o T7 /(y-( s )) ds ]) (/eC+[0,l]), (5.46) 

where Tl is the invariant law of (y(i))t>o from Corollary 15.41 and t 7 is an exponential time 
with mean 7, independent of (y(t))t>o- Setting g(x, y) := —f(y) in (|5.45|) . using the ergodicity 
of (y(t))t>o (see Corollary 15. 4j) . we find that for each z G [0, 1] and t > 0, 

n(dywr e -J"o*/(yW)d8i = lim /" p^ y(r ) € dyi^r e -/o*^y( a )) ds i 

r-+°oJ ■ (5.47) 
= lim S2S?l(x,z). 

r— >oo 

It follows from Lemma 15.131 that for each fixed r,t, and 0, the function x 1— » S^iSf l(x, z) is 
nondecr easing if / is nonincr easing, and nondecreasing and convex if / is nonincr easing and 
concave. Therefore, taking the expectation over the randomness of r 7 , the claims follow from 
(ETi6l and (ICT7I) . I 

We still need to prove Lemmas 15. 121 and 15. 131 

Proof of Lemma 15. 121 It is easy to see that the operator B from Q5.43JI is densely defined, 
satisfies the positive maximum principle, and maps the constant function 1 into 0. Therefore, 
by Hille-Yosida (|5.41|) . we must show that the range 1Z(1 — aB) is dense in C([0, l] 2 ) for some, 
and hence for all a > 0. Let V n denote the space of polynomials on [0, l] 2 of n-th and lower 
order, i.e., the space of functions / : [0, l] 2 — » M of the form 

f(x, y) = a ki x k y l with a k j = for k + I > n. (5.48) 

fcj>0 
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Set Too '■= Un'Pn- It is easy to see that B maps the space V n into itself, for each n > 0. Since 
each V n is finite-dimensional, a simple argument (see |EK86| Proposition 1.3.5]) shows that 
the image of Voo under 1 — aB is dense in C([0, l] 2 ) for all but countably many, and hence for 
all a > 0. I 

As a first step towards proving Lemma 15.131 we prove: 

Lemma 5.14 (Smooth solutions to Laplace equation) Let a > 0, g G C^([0, 1]), g < 0, 
v G C([0, l] 2 ), a^d assume that u € C^°°^([0, l] 2 ) solves the Laplace equation 

{l-a{B + g))u = v. (5.49) 

(a) If g € Ci + , f/ien v G C+ H Ci+ implies u G C+ n Ci+. 

(b) If g £ C\ + n C2+, i/ien u G C + n C\ + n C2+ implies u G C+ D Ci+ D C2+. 

Proof Let := ^u, tt 3 ^ := 

dxdy u i e ^ C- denote the partial derivatives of u and similarly for 
v and whenever they exist. Set c := ^. Define linear operators 5' and B" on C([0, l] 2 ) with 
domains V(B') = V(B") := C(°°)([0, l] 2 ) by 

B':=y(l-y)^ + {c(x-y) + 2(\-y))l r 

B" := y(l -y)$ + (c(x - y) + 4(1 - y)) 

Then 

JLBu = (B'-c) U y, 9_B'u=(B" -c-2) U y, 
J^Bu = Bu x + cuy, ^B'u = B'u x + cuy. 



(5.50) 



(5.51) 



Therefore, it is easy to see that 

(i) (l-a(B' - c + g))u y =v y + ag y u, 

(ii) (1 — a(B + g))u x = v x + a(cu y + g^u), 

(iii) (1 - a{B" - 2c - 2 + & ))tiW = u TO + a(2gy U y + 5 ra u), (5.52) 

(iv) (1 - a(B' - c + g))v?v = v x y + a(cuW + + + ftp), 

(v) (1 - a(B + <?)K X = ^ x + a(2cu xy + 2g x u x + 5 xx u), 

where in (i) and (ii) we assume that v G C^QO, l] 2 ) and in (iii)-(v) we assume that v G 
C( 2 )([0, l] 2 ). By Lemma 15.121 the closure of the operator B generates a Feller processes 
m [0, l] 2 . Exactly the same proof shows that B' and B" also generate Feller processes on 
[0, l] 2 . Therefore, by Feynman-Kac, u is nonnegative if v is nonnegative and u y ,...,u xx 
are nonnegative if the right-hand sides of the equations (i)-(v) are well-defined and non- 
negative. (Instead of using Feynman-Kac, this follows more elementarily from the fact that 
B,B', and B" satisfy the positive maximum principle.) In particular, if g v ,g x > and 
v G CW([0, l] 2 ) 5 v,v y ,v x > 0, then it follows that u,u y ,u x > 0. If moreover g yy ,g xy ,g xx > 
and v G ([0, l] 2 ), > 0, then also u vv ,u xy ,u vv > 0. I 

In order to prove Lemma 15.131 based on Lemma 15.141 we will show that the Laplace equation 
(|5.49|) has smooth solutions u for sufficiently many functions v. Here 'suffiently many' will 
mean dense in the topology of uniform convergence of functions and their derivatives up to 
second order. To this aim, we make C^ 2 ^([0, l] 2 ) into a Banach space by equipping it with the 
norm 

IMI(2) := IN + Kll + ll« x ll + \\ uVV \\ + 2 IK1 + ( 5 - 53 ) 

Here, to reduce notation, we denote the supremumnorm by ||/|| := ||/||oo- Note the factor 2 
in the second term from the right in (|5.53|) . which is crucial for the next key lemma. 
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Lemma 5.15 (Semigroup on twice diffferentiable functions) The closure inC^([0, l] 2 ) 
of the operator B generates a strongly continuous contraction semigroup on C^ 2 \[0, l] 2 ). 

Proof We must check the conditions (i)-(iii) from (|5.39|) . It is well-known (see for example 
EK86, Proposition 7.1 from the appendix]) that the space Voo of polynomials is dense in 
C( 2 )([0, l] 2 ). Therefore V(B) = C(°°)([0, l] 2 ) is dense, and copying the proof of Lemma EH 
we see that 72. (1 — aB) is dense for all but countably many a. To complete the proof, we must 
show that B is dissipative, i.e., that 

||(1 - eB)u\\ m > \\u\\ (2) (e>0,ue C<-°°\[0, l] 2 )). (5.54) 

Using (|5.51jl . we calculate 



U 1 


- eB)u = 


(1 


- e(B' - 


c))u\ 




-eB)u = 


(1 


- eB)u x 


— ECU y , 


w { 


- eB)u = 


(1 


- e(B" - 


2c-2))u yy , 


dxdy \ 


- eB)u = 


(1 


- e{B' - 


c))u xy — ecu 




-eB)u = 


(1 


- eB)u xx 


' - 2Ecu xy . 



(5.55) 



Using the disipativity of B, B' , and B" with respect to the supremumnorm (which follows from 
the positive maximum principle) we see that ||(1 — e(B' — c))u y \\ = (1 +ec)||(l — j^.B)u y \\ > 
(1 + ec) 1 1 it 37 1 1 etc. We conclude therefore from Q5.55JI that 

||(1 - £B)u\\ {2) > ||(1 - eB)u\\ + ||(1 - e(B' - c))u y \\ + ||(1 - eB)u x \\ - Ec\\u y \\ 

+ 11(1 - e{B" - 2c - 2))u yy \\ + 2||(1 - e(B' - c))u xy \\ - 2Ec\\u yy \\ 
+||(l- £j B)n-||-2 £C ||^|| 

5.5o 

> ||u|| + (1 +£c)||^|| + ||u x || -Ec\\u y \\ 
+{l + E{2c + 2))\\u yy \\+2{l+Ec)\\u xy \\-2Ec\\u yy \\ 
+ \\u xx \\ -2ec\\u xy || > ||«|| (2) 

for each e > 0, which shows that B is dissipative with respect to the norm || • ||( 2 ). ■ 

Proof of Lemma EU Let g G C^([0, l] 2 ). Then u i — ^ gu is a bounded operator on both 
C([0, l] 2 ) and C (2) ([0, l] 2 ), so we can choose a A > such that 

\\gu\\ < X\\u\\ and ||<?ii||(2) < ^IMI(2) (5.57) 

for all u in C([0, l] 2 ) and C (2) ([0, l] 2 ), respectively. Put g := g - A. By Lemma I5T21 B + g 
generates a strongly continuous contraction semigroup (<Sf)t>o = (e - Sf)t>0 on C([0, l] 2 ). 
Note that 1Z(1 — a(B + g)) is the space of all v G C([0, l] 2 ) for which the Laplace equation 
(1 — a(B + g))u = v has a solution u G C^°°\[0, l] 2 ). Therefore, by Lemma l5.14( for each 
a > 0: 

(i) If g G Ci+, then (1 - a(B + g))' 1 maps 7£(1 - a(B + g)) n C+ n C 1+ into C+ n Ci+. 

(ii) If 5 G Ci+ n C 2+ , then (1 - a{B + g))~ l maps - a(£ + g)) flC+n Ci + n C 2+ 
into c + n Ci+ n C 2+ ■ 

(5.58) 
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By Lemma 15.151 the restriction of the semigroup (Sf)t>o to C^([0, l] 2 ) is strongly continuous 
and contractive in the norm || • |b 2 ). Therefore, by Hille-Yosida (|5.39|) . 1Z(1 — a(B + g)) is 
dense in C (2) ([0, l] 2 ) for each a > 0. It follows that 11(1 - a(B + gj) D C + D Ci + is dense in 
C + H C\ + and likewise 1Z(1 — a(B + g)) n C + n C% + n C 2 + is dense in C + n C\ + n C 2+ , both in the 
norm || • ||( 2 ). Note that we need density in the norm || • ||( 2 ) here: if we would only know that 
1Z(1 — a(B + g)) is a dense subset of C([0, l] 2 ) in the norm || • ||, then 1Z(1 — a(B + g)) nC + nCi + 
might be empty. By approximation in the norm || • ||( 2 ) it follows from H5.58j) that: 



(i) If g G Ci + , then (1 — a(B + g)) 1 maps C + D C\ + into itself. 

(ii) If g G C 1+ n C 2+ , then (1 - a(B + g))~ l maps C + n C 1+ n C 2+ into itself. 

Using also continuity in the norm || • || we find that: 

(i) If g G Ci + , then (1 — a(B + g))~ l maps C + fl Ci + into itself. 



(ii) If g G Ci + D C 2 +, then (1 — a(B + 5)) 1 maps C + n Ci + n C 2+ into itself. 
For e > let 



n! 

n=0 



and therefore, by (|5.6U|) . for each t > 0: 

(i) If 5 G Ci+, then e Get maps C+ Pi Ci+ into itself. 

(ii) If g G C\ + n C 2+ , then e Get maps C + n Ci+ n C 2+ into itself. 

Finally 

e~ xt S? u = S?u = 



so ()5.63|) implies that for each t > 0: 



(i) If g G Ci + , then Sf maps C + n C\ + into itself. 

(ii) If g G Ci+ n C 2+ , then Sf maps C + n C\ + D C 2+ into itself. 



(5.59) 



(5.60) 



G £ := £ - 1 ((l-e( J B + 5))- 1 -l) (5.61) 
be the Yosida approximation to B + g. Then 



H J2 -,{l-e(B + ~g)r n (t>0), (5.62) 



(5.63) 



e - At 5f u = 5fu = lim e^u (t > 0, u G C([0, l] 2 )), (5.64) 



(5.65) 



Using the continuity of Sf in g (which follows from Feynman-Kac (|5.42jl ) we arrive at the 
statements in Lemma 15.131 I 



6 Convergence to a time-homogeneous process 
6.1 Convergence of certain Markov chains 

Section is devoted to the proof of Theorem 13.21 In the present subsection, we start by 
formulating a theorem about the convergence of certain Markov chains to continuous-time 
processes. In Section 16.21 we specialize to Poisson-cluster branching processes and superpro- 
cesses. In Section l6~3l finally, we carry out the necessary calculations for the specific processes 
from Theorem 13. 21 
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Let E be a compact metrizable space. Let A : V{A) — > C(E) be an operator defined on a 
domain T>(A) C C(E). We say that a process y = (yt)t>o solves the martingale problem for 
A if y has sample paths in T>e[0, oo) and for each / £ T)(A), the process (M/ ) t >o given by 

M/ := /(y t ) - /* A/(y a )ds (t > 0) (6.1) 
J o 

is a martingale with respect to the filtration generated by y. We say that existence (unique- 
ness) holds for the martingale problem for A if for each probability measure \x on E there is 
at least one (at most one (in law)) solution y to the martingale problem for A with initial 
law £(yo) = V- If both existence and uniqueness hold we say that the martingale problem 
is well-posed. For each n > 0, let X^> = (xf^\ . . . , X^ n ^) (with 1 < m(n) < oo) be a 
(time-inhomogeneous) Markov process in E with fc-th step transition probabilities 

P k (x,dy) = P[X^ ] £&y\x£} 1 =x] (l<fc<m(n)). (6.2) 

We assume that the are continuous probability kernels on E. Let {s^)\<k<m{n) be positive 
constants. Set 

A^f(x) := (e^r 1 [j E P k (x,dy)f(y) - f(x)) (1 < k < m(n), / G C(E)). (6.3) 



(l<*<m(n)), (6.4) 



Define := and 

/(") ._ ^ 

and put 

fe(™)(t) := max {/c : < k < m(n), 4" } < *} (* > °)- ( 6 - 5 ) 
Define processes y( n ) = (y[ n ^)t>o with sample paths in T>e[0, oo) by 

y t W :=^a (t ) («>0). (6-6) 

By definition, a space A of real functions is called an algebra if A is a linear space and f,g £ A 
implies fg G A 

Theorem 6.1 (Convergence of Markov chains) Assume that C(X^ n) ) => \i as n — > oo 

for some probability law fi on E. Suppose that there exists at most one (in law) solution to the 
martingale problem for A with initial law fi. Assume that the linear span ofV(A) contains an 
algebra that separates points. Assume that 

m(n) 

(i) lim V 4 n) = oo, (ii) lim sup e£ n) = 0, (6.7) 

k=l k: t["><T 

and 

lim sup \\A^ ) f-Af\\ oo = (fe-D(A)) (6.8) 

n ^°° k: t ™<T 

for each T > 0. Then there exists a unique solution y to the martingale problem for A with 
initial law fi and moreover £(y( n )) =4> C(y), where =>- denotes weak convergence of probability 
measures on T>e[0, oo). 
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Proof We apply |EK861 Corollary 4.8.15]. Fix / G V(A). We start by observing that 



f(4 U) ) - E ^ } A (n) /(^|) (0 < k < m{n)) 



.(«) /t( n ) ft v( n ) 



(6.9) 



i=l 



is a martingale with respect to the filtration generated by and therefore, 

fc(«)(t) 

/(yl n) ) -E ^A w /(yS ) (* > o) 
• 1 »-i 

is a martingale with respect to the filtration generated by y ( n ) . Put 



and set 



and 



» 



m(n) J 



d n) :=/(yj n) )+ f ^ds (i>0). 



'lt\ 

Then we can rewrite the martingale in (|6.10|) as 

rt 



i n) ds. 



(6.10) 

(6.11) 
(6.12) 

(6.13) 
(6.14) 



o 



By |EK86| Corollary 4.8.15] and the compactness of the state space, it suffices to check the 
following conditions on (j)^ and £( n ); 

(i) sup sup£[|£ t (n) |] < oo, 

n>N t<T 

(ii) sup sup£[|4 n) |] < oo, 

n>N t<T 



(iii) lim E 

n— >oo 



(4 n) -/(y? 3 ))II^(y& 



» 



n) 1 



0. 



i=l 



(iv) lim E (4 n) - Af(y^)) J] h t (y^ 



(n)^ 



i=l 



(6.15) 



(v) lim 

n— >oo 



S up |d n) -/(yi n) )|l =o, 

■teQn[o,T] J 
(vi) sup ^[H^^^llp^] < oo for some p £ (1, oo], 



n>N 



for some N > and for each T > 0, r > 1, < si < ■ ■ ■ < s r < T, and h\, . . . , h r € H C C(E). 
Here 7i is separating, i.e., J /id^, = J hdv for all h £ TC implies \x = v whenever fj,, v are 
probability measures on E. In (vi): 



\9\\p,T 



T \Vp 



n 



(1 < p < oo) 



(6.16) 



and ||g||oo,T denotes the essential supremum of g over [0, T]. 
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The conditions (|6.15() (i)-(vi) are implied by the stronger conditions 

(i) Jim snp |kJ n) -/(yJ ri) )|L = 0, 



— 0<t<T ,r~ -yrj\ 

(ii) lim sup ||4 n) -A/(y i (n) )|| oo = 0, 



where we denote the essential sup re mum norm of a real- valued random variable X by || JTUqq : — 
mf{K > : \X\ < K a.s.}. Condition (IQ71) (ii) is implied by (jHHJ) (i) and ffTBjl . To see that 
also (|6.17() (i) holds, set 

M n := sup H^IL, (6.18) 

0<t<T 

and estimate 

sup \\^ n) - /(y! n) )L < M n sup{4 n) : 1 < k < m(n), < T}. (6.19) 

0<t<T 

Condition (|6. 17|) (ii) implies that lim sup n M n < oo and therefore the right-hand side of ()6.19|) 
tends to zero by assumption (|6,7[) (ii). I 



6.2 Convergence of certain branching processes 

In this section we apply Theorem 16. II to certain branching processes and superprocesses. 

Throughout this section, E is a compact metrizable space and A : T>(A) —* C{E) is a linear 
operator on C(E) such that the closure A of A generates a Feller process ^ = (&)t>o in E with 
Feller semigroup (P t ) t >o given by P t f(x) := E x [f(^)} (t > 0, / G C{E)). 

Let a G C+(E) and /?, / G C(E). By definition, a function 1 1— > « t from [0, oo) into C(-E) is 
a classical solution to the semilinear Cauchy problem 

| §- t u t = Au t + (3u t -au 2 t {t>0), 
{ u = f 

if t ut is continuously differentiable (in C(E)), ut G *D(A) for all t > 0, and (|6,2()|) holds. 
We say that u is a mild solution to (|6,2U[) if 1 1— ► u t is continuous and 

u t = P t f+ [ P t - s (pu s -au 2 s )ds (t>0). (6.21) 

JO 



Lemma 6.2 (Mild and classical solutions) Equation \6.2U\) has a unique Cj r {E)-valued 
mild solution u for each f G C+(E), and f > implies that ut > /or all t > 0. If moreover 
f G *D{A) then u is a classical solution. For each t > 0, ut depends continuously on f G C+(i£). 

Proof It follows from |Paz831 Theorems 6.1.2, 6.1.4, and 6.1.5] that for each / G C(£), (l6~277|) 
has a unique solution («t)o<t<T up to an explosion time T, and that this is a classical solution 
if / G ^(^4)- Moreover, iif depends continuously on /. Using comparison arguments based 
on the fact that A satisfies the positive maximum principle (which follows from Hille-Yosida 
(|5.41[l ) one easily proves the other statements; compare |FS04[ Lemmas 23 and 24]. I 

We denote the (mild or classical) solution of (|6.20|) by Utf := Uf, then Ut : C+(E) — > C+(E) 
are continuous operators and U = {Ut)t>o is a (nonlinear) semigroup on C + (E). 
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Since E is compact, the spaces {fj, G M.(E) : n(E) < M} are compact for each M > 0. In 
particular, M.{E) is locally compact. We denote its one-point compactification by j^i(E) 00 = 
M(E) U {oo}. We define functions Fj G C{M{E) 00 ) by F/(oo) := and 

F f {p) := e-^' /> (/ G C + (F), / > 0, // G M(F)). (6.22) 
We introduce an operator £/ with domain 

25(0) := {-F/ : / G D(A), / > 0}, (6.23) 

given by QFf(oo) := and 

gF f ( fJ ,):=-(v,Af + (3f-af)e-(t J '>f) (^GM(E)). (6.24) 
Note that QF f G C{M.{E) 00 ) for all F/ G £>(£). 

Proposition 6.3 ((A, a, /?)-superprocesses) TTie martingale problem for the operator Q is 
well-posed. The solutions to this martingale problem define a Feller process y = (3^)*>0 in 
«M(F)oo with continuous sample paths, called the (A, a, (3)-superprocess. If y$ = oo then 
y t = oo for all t>0. Ify = fj,£ M(E) then 

£M[ e -Cyt,/}] = e -(n,U t f) (/gC+(F)). (6.25) 

Proof Results of this type are well-known, see for example |KK86I Theorem 9.4.3], |h'it,88j . 
and |EK9U Theoreme 7]. Since, however, it is not completely straightforward to derive the 
proposition above from these references, we give a concise autonomous proof of most of our 
statements. Only for the continuity of sample paths we refer the reader to |Fit881 Corol- 
lary (4.7)] or ER91, Corollaire 9]. 

We are going to extend Q to an operator Q that is linear and satisfies the conditions of 
the Hille-Yosida Theorem Q5.41J1 . For any 7 G C + (E) and \i G A4(E), let Clust 7 (/x) denote a 
random measure such that on {7 = 0}, Clust 7 (/i) is equal to and on {7 > 0}, Clust 7 (/i) is a 
Poisson cluster measure with intensity ^/x and cluster mechanism Q(x, •) = C(tj^6 x ), where 
is exponentially distributed with mean 7(3;). It is not hard to see that 

£[ e -(Clust 7 (//),/}] = e -(n,Vjf) (/ g C(E), f > 0), (6.26) 

where V 7 /(x) := (j^y + 7(^)) _1 - Note that since V 7 1 is bounded, the previously mentioned 
Poisson cluster measure mentioned above is well-defined. By definition, we put Clust 7 (oo) := 
00. 

Define a linear operator Q a on C(M.(E)) OQ ) by 

G a F(fx) := lime- 1 (E[F(OxiBt ea (jj.))] - F(jj)) (6.27) 

with as domain T>(Q a ) the space of all F G C(M(E) OQ ) for which the limit exists. Define a 
linear operator Qp by 

^F( M ) := lim e- 1 (F((l + e/3) M ) - Fiji)) (6.28) 
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with domain V{Gp) := C(7W(^)) 00 ). Define P t * : M^E)^ -» MiE)^ by (P t >,/) := <//,i> t /> 
(t>0, f eC(E), fie M(E)) and P t *oo := oo (t > 0). Finally, let ^ be the linear operator 
on C(.M(-E))oo) defined by 

QjFip) := lime-^F^V) - (6.29) 
with as domain T>(Gj) the space of all F for which the limit exists. Define an operator Q by 

g ■= g a + g p + g^, (6.30) 

with domain V(Q) := T>(G a ) n If / G T>(A), f > 0, and F/ is as in (f6~2"2l. then it is 

not hard to see that GFf (oo) = and 

GFf{n):=-{^Af + {3f-af 2 )e-^f) (f,eM(E)). (6.31) 

In particular, extends the operator £/ from (|6.24|) . Since 2?(.A) is dense in C(E), it is easy to 
see that {F f : f G P(3), / > 0} is dense in C(A^(F) 00 ). Hence X>(<?) is dense. Using (I7T2TJ) - 
1)6.29(1 it is not hard to show that G satisfies the positive maximum principle. Moreover, by 
Lemma IO for / G V(A) with / > 0, the function t i-> F Utf from [0, oo) into C(M(E) 00 ) is 
continuously differ entiable, satisfies Fy t f 6 T>(G) for all t > 0, and 

|i^ t/ = (t > 0). (6.32) 

From this it is not hard to see that G also satisfies condition (|5.41|) (ii), so the closure of 
G generates a Feller semigroup (St)t>o on C(M(E) 00 ). It is easy to see that SfFf = Fy t f 
(t > 0). By EKSl Theorem 4.2.7], this semigroup corresponds to a Feller process y with 
cadlag sample paths in M.(E)oo. This means that E fJ, [Ff(yt)] = Fu t f(fj) for all / G V(A) 
with / > 0. If jj, = oo this shows that 3^ = oo for all t > 0. If /i G M(E) we obtain (|6.25l) for 
/ G P(-A), / > 0; the general case follows by approximation. I 

Now let (q £ )e>o be continuous weight functions and let (Q e ) e >o be continuous cluster mecha- 
nisms on E. Assume that 

Ze(x) := f Q £ (x,d X )(xA) <oo (x G E) (6.33) 

and define probability kernels K £ on E by 

J K £ (x, dy)f(y) := -^rr J Q E (x, d X )( X , f) (f G B{E)). (6.34) 

For each n > 0, let (£j, n ^)i<fe< m ( n ) (with 1 < m(n) < oo) be positive constants. Let = 
(Xn , ■ ■ ■ , X^ n ) .) he & Poisson-cluster branching process with weight functions q („),..., q („) 

and cluster mechanisms Q („) , . . . , Q ( n ) . Define tj^ and k^ n \t) as in (|6.4j) - (|6.5j) . Define 
processes by 

# : = (*>0). (6-35) 
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Theorem 6.4 (Convergence of Poisson-cluster branching processes) Assume that 
C(Xq ) p as n — > oo for some probability law p on Ai(E). Suppose that the constants 
ejj. fulfill j6.7\) . Assume that 

(i) g e (x) J Q £ (x,d x )(xA) = 1 + e(3(x) + o(e), 

(ii) g £ (x)y Q £ (x,d X )(x,l) 2 =^2a(x) + (e), (6.36) 

(iii) Qe(x) / Q e (x,dx)(x, l) 2 l{( x ,i>>«5} =o(e) 



/or each 5 > 0, and 

J K £ (x, dy)f{y) = f(x) + eAf(x) + o(e) (6.37) 

for each f G TJ(A), uniformly in x as e — > 0. Then £(y( n ') C(y), where y is the (A, a, j3)- 
superprocess with initial law p. 

Here =>■ denotes weak convergence of probability measures on T>j^fm[0, oo). 

Proof We apply Theorem 16.11 to the operator Q, where we use the fact that if we view 
Mi(T>m(e) [0) °°)) as a subspace of A^i^^i^^ [0, oo)) (note the compactification) , equipped 
with the topology of weak convergence, then the induced topology on Mi(Vj^^[0, °°)) is 
again the topology of weak convergence. 

By Proposition 16.31 solutions to the martingale problem for Q are unique. Since FfF g = 
Ff +g and T>(A) is a linear space, the linear span of the domain of Q is an algebra. Using the 
fact that T>(A) is dense in C(E) we see that this algebra separates points. Therefore, we are 
left with the task to check (|6.8|) . 

Define U £ : C + {E) C + {E) by 

U £ f(x) := q £ (x) J Q £ (x,d X ){l- e-(xJ)) ( x G E , feC+[0,l], f > 0, s > 0), (6.38) 
and define transition probabilities P £ (p,,dv) on _M(£)oo by P £ (oo, ■) := 6^ and 

P e (ji,dj/)e-( u >f) = e-(v,UJ). (6.39) 



/ 



We will show that 



hm || e - 1 (^/-/)-(^/ + /3/-a/ 2 )|| oo = (feV(A), f > 0). (6.40) 

Together with (|6.39|) this implies that 

J P £ (p,dv)F f (v) = Ff (p)+egF f (p) + o(e) (/ G / > 0), (6.41) 

uniformly in € M(E) QO as e — > 0. Therefore, the result follows from Theorem 16. II 
It remains to prove (|6.4U|) . Set g{z) := 1 — z + ^z 2 — e~ z (z > 0) and write 

U £ f{x)=q £ {x) [ Q £ (x,dx){(x,f)-h(xJ) 2 +g((x,f))). (6.42) 
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Since 

rv r x 



g (z) = dy dx dte' 1 (z > 0), (6.43) 
Jo Jo Jo 



it is easy to see that g is nondecreasing on [0,oo) and (since < e * < 1 and f Q x dte 1 < 1) 



2 

Using these facts and (|6.36[) (ii) and (iii), we find that 



< g(z) < \z 2 A \z z (z > 0). (6.44) 



Qe{x) J Qe(x,dx)g((x,f)) 

< ll/lloo? e (x){ I Q E (x,dx)g((xA))l{( x ,i)<5} + I Qe(x,d x )g((xA))l{( x ,i)>8}} 

< ll/ll«o*(*){»*^ 

= \5\\f\\ 00 {e2a{x) + o(e))+o{e). 

(6.45) 

Since this holds for any 5 > 0, we conclude that 

q £ (x) J Qe(x,dx)g((x,f}) = o(e) (6.46) 

uniformly in x as e — > 0. By (|6.36[) (i) and (j6,37|) . 

qe(x) j Qe(x,d X )(xJ) = (qs(x) J Q £ (x, d X ) (x, 1>) ( J K £ (x , dy) f (y)) 

= (1 + e(3(x) + o(e)) (f(x) + eAf(x) + (e)) ' ' (6-47) 

= f(x) + eP(x)f(x) + eAf(x) + o(e). 

Finally, write 

qe(x) [ Qe(x,d X )(x,f) 2 

J , (6.48) 

= Qe(x) J Q £ (x,dx){(xJ(x)) 2 + 2(x,f(x))(x,f-f(x)) + (x,f-f(x)) 2 ). 

Then, by (IQHll (ii), 

q £ (x) [ Q £ (x,d X )(x,f(x)) 2 = f(x) 2 {e2a(x) + o(e)). (6.49) 



< (? e (x) y Q e (x,d X )(x,/-/(x)) 2 ) / Q e (x,dx)(x,/W) : 

< (o{e){2a{x)e + o(e))) 1/2 = o(e). 



We will prove that 

? e (ar) / Qe(x, dx)(x, / - /(x)> 2 = o(e). (6.50) 



Then, by Holder's inequality, (|6.36|) (ii), and (|6.5Uj) . 

\q £ (x) J Q £ (x,d X )(xJ-f(x))(x,f(x))\ 

.1/2/ _ /■„ , 2 y/2 (6 51) 
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Inserting ((OS)) . (KPtUI) and (I63T1) into we find that 

qe(x) ! Qe(x,dx)(x,f) 2 = e2a(x)f(x) 2 +o(e). (6.52) 



Inserting (HOIJ), (|6~4T)) and (jB~52l into (|6~4^|) . we arrive at (|6~iT)|) . We still need to prove 
(|6.50j) . To this aim, we estimate, using (|6.47|) . 



(6.54) 



qe(x) J Qe(x,dx){x,f ~ f(x)) 2 l{( x ,i)<8} 

<5||/-/(x)||oog £ (x) J Qe(x,dx)(X,f-f(x)) ^ 

= 5\\f-f(x)\\ 00 {eAf(x) + o(e)) 

and, using ()6.36j) (hi), 

Qs(x) J Qs(x,dx)(x,f ~ f(x)) 2 l{( x ,i)>8} 

< ||/- f(x)\\ooq E (x) J Q E (x,dx)(x,l) 2 l{( Xl i»5} = o(e). 

It follows that 

g e (z) y Q e (x,dx)(x,/-/(x)) 2 <fe||/-/(x)|| 00 ^/(x) + (e) (6.55) 
for any 5 > 0. This implies (|6.50j) and completes the proof of (|6.40|) . I 

6.3 Application to the renormalization branching process 
Proof of Theorem 13.21 (a) For any /q, . . . , fk £ C + [0, 1] one has 

E[e~ (X-n, fo) . . . e — (X- n +k, fk/] 

= E[e~( X - n ' f°) ■■■ e-^-n+k-Xifk-x + U ln _ k f k )^ 
= ...= E[e~^ X - n ^k}^ 

where we define inductively 

g :=fk and g m+ i := fk-m-i + U ln _ k+m g m . (6.57) 

By the compactness of [0, 1] and Corollary 15.101 the map (7, /) 1— > U^f from (0, 00) x C+[0, 1] 
to C+[0, 1] (equipped with the supremumnorm) is continuous. Using this fact and (|6.56|) we 
find that 

E\e~( X - n > ■ ■ ■ e _ ^- n + fc ' f k h > ^^"(^-n, /o) . . . e _ (3 ; _ n+ fc> A)i (6.58) 

^ J n— >oo '- 

Since fi,- ■ ■ , fk ar e arbitrary, (|3.12|) follows. I 

Proof of Theorem 13.21 (b) We apply Theorem 16.41 to the weight functions g 7 and cluster 
mechanisms Q 7 from Q3.8JI and to j4wf = ^(1 — x )~§^? with domain T>(A^p) = C^[0, 1], and 
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a = j3 = 1. It is well-known that Awf generates a Feller semigroup EK86, Theorem 8.2.8]. 
We observe that 

J Q,(xAx){xJ) = E[2 £ f{yl{-t))]=2E[r 1 ]E[f(yim] = 7 j TZ(dy)f(y), (6.59) 

where Y% is the equilibrium law of the process y2 from Corollary 15.41 It follows from ([5.24ft 
that 

(i) / ri(dy)(y-x) = 0, 



2 7X(1 - x) 



(ii) jTZ(dy)(y-x) j + . 

(iii) fvi(dy){y- X f = 0{ 1 % 



(6.60) 



uniformly in x as 7 — » 0. Therefore, for any 5 > 0, 

(i) / n(dy)(i/-a;) = 0, 

r 2( d y)(y - 2;) 2 = 7^(! - x ) +0(7), 
(iii) y T2(dy)l {ly _ xl>5} = o('y), 

uniformly in x as 7 — > 0. Consequently, a Taylor expansion of / around x yields 

'T2(dy)f(x) = f{x) + 1 \ x {l-x)£ I f{x) + 0(7) (/ e C( 2 )[0,1]), 



(6.61) 



(6.62) 



uniformly in x as 7 — > 0. (For details, in particular the uniformity in x, see for example 
Proposition Swa99] B.l.l].) This shows that condition (|6.37j) is satisfied. Moreover, 

J Qj(x,d X )(xA) =E[2 Tj ] = 7 , 

/f'OO 
Q 1 (x,dx)(xA) 2 = E[(2r 7 ) 2 } = j z 2 \e- z ^dz = 2 7 2 , (6.63) 

/9*OG 
Q y (x,d X )(xA) 3 = E[(2r J f] = j z^e^dz = 6 7 3 , 

which, using the fact that g 7 = (- + 1), gives 

q-y J Q 7 (x,dx){x, 1 ) = 1 + 7, 

<7 7 y Q 7 (x, dx)(x. I) 2 = 27 + 0(7), (6.64) 
q y J Q y (x,d X )(xAf =0(7). 
This shows that (|6.36j) is fulfilled. In particular, 

q~ ( J Qj(x, d X )(x, 1) 2 1{( X ,1»<5} < *~ X <h J d x)(x, I) 3 = o(7) (6.65) 

for all 5 > 0. I 
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7 Embedded particle systems 

In this section we use embedded particle systems to prove Proposition 13.51 An essential 
ingredient in the proofs is Proposition 17.151 (a), which will be proved in the Section |HJ 

7.1 Weighting and Poissonization 

Proof of Proposition 13.31 Obviously q% G C + (E h ) for each k = 1, . . . ,n. Since h G C + (E) 
and h is bounded, it is easy to see that the map \i \—* h\x from A4(E) into A4(E h ) is continuous, 
and therefore the cluster mechanisms defined in (|3.21|) are continuous. Since 

U " f{x) = f^H 1 - e ~ (hZxJ) ] = m h$ X) ( x ^ E h , f G B + (E h )), (7.1) 

formula (|3.22|) holds on E h . To see that (j3.22|) holds on E\E h , note that by assumption 
Ukh < Kh for some K < oo, so if x G E\E h , then Ll k h(x) = 0. By monotonicity also 
Uk{hf){x) = 0, while hU k f(x) = by definition. Since snp xeE hU k l(x) = sup xeE h < 

K < oo, the log-Laplace operators satisfy (|3,3|) . If A' is started in an initial state A?o, 
then the Poisson-cluster branching process X h with log-Laplace operators . . . started 
in Xq = hX$ satisfies 

E r e -(hX k ,fh =E [ e -(X ,Ui ° • • • °Z4(/i/))j 

= E[e -(X ,hU?o...oUj;(f))] =E[e-( X tf)] (f£B + (E h )), ^ 

which proves (|3.23j) . I 
Proof of Proposition 13.41 We start by noting that by (|3.2j) . 

U k f(x) = q(x)E[l - e-^J)] = qk ( x )P[Pois(fZ%) + 0] (x € E, f G £+(£)). (7.3) 
Into (|3.24|) . we insert 

p[PdB(fc2*)e.] 

= P[Pois(fcZ*) G • | Pois(hZ k ) / 0] PpPois(/»Z*) / 0] + 5 P[Pois(hZ%) =0]. 

Here and in similar formulas below, if in a conditional probability the symbol Pois( • ) occurs 
twice with the same argument, then it always refers to the same random variable (and not to 
independent Poisson point measures with the same intensity, for example). Using moreover 
(|7.3j) we can rewrite (|3,24j) as 

Q h k (x, • ) = ^P[Pois(^) G • | Pois(^) + 0] + HX) ~^ H{X) Soi ■ )• (7-5) 

In particular, since we are assuming that h is Wfc-subharmonic, this shows that Q k (x, • ) is 
a probability measure. Let X h be the branching particle system with offspring mechanisms 
Qi, . . . ,Q k . Let Zx ,k be random variables such that C{Zx' k ) = Q k (x, •). Then, by Q3.18JI . 
lt3~54l . dim , and Q, 



UZf{x) = P[Thm f (Z^ k ) * 0] = |^P[Thin / (Pois(/ i ^)) * 0] 

= T$ P[Fois{hfz * ) * 0] = l^x) Uk{hf){x) {x e Eh) - 



(7.6) 
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If x G E\E h , then U k {hf)(x) < U k (h)(x) < h{x) = =: W h {f){x). This proves (l3~23j) . To 
see that Q\ is a continuous offspring mechanism, by |Kal761 Theorem 4.2] it suffices to show 

that J Q^(x, dv)e~ & 3) i s continuous for all bounded g G C + {E h ). Indeed, setting / := 
l-e-9, one has J Q^(x,du)e~^9) = J Q h k { x ,dv)(l- f f = 1-U%f(x) = l-U k (hf){x)/h(x) 
which is continuous on E h by the continuity of q k and Q k . 

To see that also (|3~26|) holds, just note that by (|3T9|) . (jH2SJ), and (|33|) . 

^ (P ° is(M) [Thi n/ (^) = 0] = P[Thm^ ... ^ / (Pois(M) = 0] 

= prPois((/»E/f o • • • o U*f)n) = 0] = P[Pois((Wi o • • • o U n (hf))ii) = 0] (7.7) 
= P M [Pois(/i/A' n ) = 0] = P /i [Thin / (Pois(/iA' n )) = 0]. 

Since this formula holds for all / G B{ 01 ](E h ), formula (|3.26j) follows. I 

Remark 7.1 (Boundedness of h) Propositions 13.31 and 13.41 generalize to the case that h 
is unbounded, except that in this case the cluster mechanism in (|3.21|) and the offspring 
mechanism in ()3.24j) need in general not be continuous. Here, in order for (|3.22|) and (|3.25|) 
to be well-defined, one needs to extend the definition of U k f to unbounded functions /, but 
this can always be done unambiguously FS03, Lemma 9]. O 

7.2 Sub- and superharmonic functions 

This section contains a number of pivotal calculations involving the log-Laplace operators U-y 
from (|3.9j) . In particular, we will prove that the functions hi i, /io,o> and /io,i from Lemmas 13.61 
13.71 and 13.81 respectively, are £/ 7 -superharmonic. 

We start with an observation that holds for general log-Laplace operators. 

Lemma 7.2 (Constant multiples) Let U be a log-Laplace operator of the form sat- 
isfying £Dj) and let f G B+(E). Then U(rf) < rUf for all r > 1, and U(rf) > rUf for all 
< r < 1. In particular, if f is U -superharmonic then rf is U -superharmonic for each r > 1, 
and if f is U -subharmonic then rf is IA -superharmonic for each < r < 1. 

Proof If X is a branching process and IA is the log-Laplace operator of the transition law from 
Xq to X\ then, using Jensen's inequality, for all r > 1, 

e -{n,U(rf)) _ j5**[ e -(<*i,r/)] = £M[( e -(-*i,/)) r ] > {E^[e~^ Xl ^)]) r = eT^ rU f) . 

(7-8) 

Since this holds for all ^ G M(E), it follows that lA{rf) < rlAf . The proof of the statements 
for < r < 1 is the same but with the inequality signs reversed. I 

We next turn our attention to the functions h\ \ and /io,o- 
Lemma 7.3 (The catalyzing function One has 

U 1 {rh hl )(x) = \^- ( 7 ,r > 0, xG [0,1]). (7.9) 

r + 7 

In particular, hi t \ is lA~ f -harmonic for each 7 > 0. 
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Proof Recall (|3.7j) - (j3,9|) . Let o~i/ r be an exponentially distributed random variable with mean 
1/r, independent of r 7 . Then 

UytrhwXx) = (i + 1)E[1 - e~ lo 1 rdt ) = (I + l)P[a 1/r < r 7 ] = (I + 1)^, (7.10) 

which yields (JUJ). I 

Lemma 7.4 (The catalyzing function ho o) One has W 7 (r/io,o) < rho,o for each 7, r > 0. 
Proof Let Tx be the invariant law from Corollary 15.41 Then, for any 7 > and / G B+[0, 1], 

Urf{x) = (I + - e" W />] < (I + /)] 

= (i + i) J B[^ T >( y 7 ( _ t /2)) d t] =(i + 7 )<r2,/) (xe[o,i\), (7 ' U) 

where we have used that r 7 is independent of yl and has mean 7. In particular, setting 
/ = rho t o and using (j5.25j) we find that £/ 7 (r/iQ,o) < r/to,o- ■ 

The aim of the remainder of this section is to derive various bounds on tl^f for / G 7io,i- We 
start with a formula for U-yf that holds for general [0, l]-valued functions /. 

Lemma 7.5 (Action of on [0, l]-valued functions) Let y2 be the stationary solution 
to \3.b)) and let r 7/ / 2 be an independent exponentially distributed random variable with mean 
7/2. Let {fli)i>i be independent exponentially distributed random variables with mean ^, in- 
dependent of yZ and t^m, and let at '■= Yli=i@i (k>0). Then 

l-U 1 f(x) = E[ H (l-/(y7(_ afe )))] ( 7>0 , /GB [Ojl] [0,l], xe [0,1]). (7.12) 

fc>0: CTfc<r 7 

Proof By Lemma 17.31 the constant function h\ ±(x) := 1 satisfies LCyhn = hi s i for all 7 > 0. 
Therefore, by Proposition 13.41 Poissonizing the Poisson-cluster branching process X with the 
density h\ \ yields a branching particle system X hl < 1 = {Xj£ , • • • , X 1,1 ) with generating 
operators U^'\ , . . . , U^' 1 , where 

ll!} 1 ' 1 f = U y f (/€B [0 ,i] [0,1], 7>0). (7.13) 

By dUIHI) and (1731) . 

[^70*0 = l-E[(l-/) Pois (^) |Pois(Zj) / 0] (/ G S [0il] [0,l], x G [0,1], 7 > 0). (7.14) 
Therefore, ()7.12j) will follow provided that 

P[Poi S (ZJ) G • I PoisOZZ) + 0] = £( J] *y2(-«r fc ))- (7-15) 

fc>0: o- fc <r 7 / 2 



Indeed, it is not hard to see that 



Pois(Zj)^ £ ^ ( _ fffc) . (7.16) 

fc>0: <J k <T 1 / 2 
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This follows from the facts that Z2 = 2 Jq i/2 <5 y 7(_ s )d,s and 

£ 5_ (Jfc ^Pois(21 ( _ r7/2i0] ). (7.17) 

k>0: a k <T 1 / 2 

Conditioning Pois(21(_ T ]) on being nonzero means conditioning on t 7 / 2 > 0*1 • Since 
T 7/2 ~~ conditioned on being nonnegative, is exponentially distributed with mean 7/2, 
using the stationarity of y2, we arrive at (|7.15|) . I 

The next lemma generalizes the duality Q5.22JI to mixed moments of the Wright-Fisher diffusion 
y at multiple times. We can interpret the left-hand side of (|7.18[) as the probability that 
mi, . . . , ra n organisms sampled from the population at times t\, . . . , t n are all of the genetic 
type I. 

Lemma 7.6 (Sampling at multiple times) Fix < t\ < ■ ■ ■ < t n = t and nonnegative 
integers mi, . . . , m n . Let y be the diffusion in \5.2U\) . Then 

n 

EV [[\yZ k } =E[y^x% (7.18) 

k=l 

where (<p s , ip s ) s &[o,t] ^ s a Markov process in N 2 started in (4>o,ipo) = (m n ,0), that jumps deter- 
ministically as 

(0s, i> 8 ) ->■ {4>s + m k , tp s ) at time t - t k (k<n), (7.19) 
and between these deterministic times jumps with rates as in \5.21\) . 



Proof Induction, with repeated application of (|5.22jl . I 

For any m > 1, we put 

h m (x) := 1 - (1 - x) m (2: €[0,1]). (7.20) 

The next lemma shows that we have particular good control on the action of on the 
functions h m . 

Lemma 7.7 (Action of on the functions h m ) Let m > 1 and let r 7 be an exponentially 
distributed random variable with mean 7. Conditional on r 7 , let (cf)' t , ip' t )t>o be a Markov process 
in N 2 , started in {(/>q,iP'q) = (m, 0) that jumps at time t as: 



t, -> (<t>'t ~ 1, *P't + 1) with rate ±<f>' t , (7.21) 



b' t , -» (# - 1, with rate - 1 

1 

7 

i>t) ~* (fit + "1, ft) With rate 1 {r 7/2 <t} 

Then the limit lim^oo tp' t =: ip'^ exists a.s., and 

Uyhmix) = £ (m ' 0) [1 - (1 - xf~] (m > 1, x e [0, 1]). (7.22) 
Proof Let y2, T7/2) an d (cfc)fe>o be as in Lemma 1731 Then, by (|7.12|) . 

U Y h m (x) = l-E[ [] {l-yl{-a k )) m ]. (7.23) 

fc>0: o-fe<T 7 / 2 



46 



Let (cfi',ip') = (4>' t ,^' t )t>o be a N 2 -valued process started in (^qjV'o) = ( m ify such that condi- 
tioned on r 7 and (cr k ) k >o, {<ft' ■,i J> ) is a Markov process that jumps deterministically as 

((ft' t ,tp' t ) -» + m,^) at time a k (k > 1 : a k < Ty/2 ) (7.24) 

and between these times jumps with rates as in (|5.21|) . Then ((ft' t ,ip' t ) — > (O,^^) as t — > oo 
a.s. for some N- valued random variable ^j^, and (|7.22|) follows from Lemma |7.61 using the 
symmetry y <-> 1 — y. Since (Tfc+i — are independent exponentially distributed random 
variables with mean one, (</>', tp') is the Markov process with jump rates as in (|7.21j) . I 

The next result is a simple application of Lemma 17.71 

Lemma 7.8 (The catalyzing function hi) The function h\(x) := x (x G [0,1]) is lA^- 
subharmonic for each 7 > 0. 

Proof Since ip'^ > 1 a.s., one has 1 - (1 - x)^'°° > x a.s. (x G [0, 1]) in (|7.22j) . In particular, 
setting m = 1 yields U 7 hi > h\. I 

We now set out to prove that hj, which is the function hoi from Lemma l3.8| is ^-super- 
harmonic. In order to do so, we will derive upper bounds on the expectation of tp'^. We derive 
two estimates: one that is good for small 7 and one that is good for large 7. 

In order to avoid tedious formal arguments, it will be convenient to recall the interpreta- 
tion of the process ((ft\ip') and Lemma 17.61 Recall from the discussion following 1)5.22(1 that 
(y2(t))tgiR describes the equilibrium frequency of genetic type / as a function of time in a 
population that is in genetic exchange with an infinite reservoir. From this population we 
sample at times — a k (k > 0, a k < T j/2) each time m individuals, and ask for the probability 
that they are not all of the genetic type II. In order to find this probability, we follow the 
ancestors of the sampled individuals back in time. Then (ft' t and ip' t are the number of ancestors 
that lived at time —t in the population and the reservoir, respectively, and E[l — (1 — x)^°°] 
is the probability that at least one ancestor is of type I. 

Lemma 7.9 (Bound for small 7) For each 7 G (0, 00) and m > 1, 

1 1 m-l 

-E^W^] < - £ =: XmiTi). (7.25) 

i=o ' 

The function Xm is concave and satisfies Xm(0) = 1 for each m > 1. 
Proof Note that 

E[\{k>0: o- k <T l/2 }\] =1 + 7. (7.26) 

We can estimate (eft', ip') from above by a process where ancestors from individuals sampled 
at different times cannot coalesce. Therefore, 

^ m ' )[^]<(l+7)S (m ' 0) [^oo], (7.27) 

where ((ft, ip) is the Markov process in (|5.21j) . Note that if ((ft, tp) is in the state (m + 1, 0), then 
the next jump is to (m, 1) with probability 

-(m+1) 1 

(7.28) 



i(m + 1) + m(m + 1) 1 + rwy 
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and to (m, 0) with one minus this probability. Therefore, 



1 + my V 1 + ?7i7 / 



+ mrj V 1 + rwy * 

: - J— (V™' ) [^oo] + l) + (l " T^— ) [^oo] (7.29) 

1 + my V / \ 1 + m~i / 

1 + m7 



By induction, it follows that 

m— 1 



£(m ' 0) w = Er^-- ( 7 - 3 °) 

i=0 ' 



Inserting this into (|7.27jl we arrive at (|7.25|) . Finally, since 

d 2 1 + 7 2i(*-l) 



57 2 l + i7 (l + »7) 3 
the function Xm is convex. 



> (i > 0, 7 > 0), (7.31) 



Lemma 7.10 (Bound for large 7) For each 7 E (0, 00) and m>l, 

m 

£ (m ' 0) [v4]<(^ + i)Efc + |- ( 7 - 32 ) 

k=l 

Proof We start by observing that ^.E[ip' t ] = ~E[(f)' t }, and therefore 

/•CO 

EWoo] = k / E[<fc]dt. (7.33) 



Unlike in the proof of the last lemma, this time we cannot fully ignore the coalescence of 
ancestors sampled at different times. In order to deal with this we use a trick: at time zero we 
introduce an extra ancestor that can only jump to the reservoir when t > r 7 and there are no 
other ancestors left in the population. We further assume that all other ancestors do not jump 
to the reservoir on their own. Let £t be one as long as this extra ancestor is in the population 
and zero otherwise, and let eft" be the number of other ancestors in the population according 
to these new rules. Then we have at a Markov process (£, <£>") started in (£q, 0q) = (1, m) that 
jumps as: 

with rate + 1)$', 
(6, <t>?) (6, <t>t + "») with rate l{r l/2 <t}, (7.34) 
(6, 0t ) (ft - 1, 0") with rate ^l{r 7/2 >t}l{^'=o}- 

It is not hard to show that (£, </>") and 0' can be coupled such that £j + 0" > 0j for all t > 0. 
We now simplify even further and ignore all coalescence between ancestors belonging to the 
process <f) that are introduced at different times. Let <p\ be the number of ancestors in the 
population that were introduced at the time o~k [k > 0). Thus, for t < o~k one has = 0, 
for t = (Tfr one has <f>^ = m, while for t > a^, the process (f>[ k ^ jumps from nton-1 with rate 
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(n + l)n. Then it is not hard to see that, for an appropriate coupling, < J2k>o-a k <r /2 'H 

(7.35) 



for all t > 0. We let £' be a process such that £g = 1 and ^ jumps to zero with rate 

1 



■1 



7 



n h 

k>0:a k <T l/2 



Then for an appropriate coupling > £ t (i > 0). Thus, we can estimate 



E[<fc]dt< / E$]dt + / S En 

^ ^ L fc>0:<, fc <r 7/2 



(AO 



dt. 



(7.36) 



Set p := inf{t > r 7 / 2 : 0^ = VA; > with < t 7 / 2 } and it := inf{i > : £,' t = 0}. Then 

f°° 3 
y B[^]dt = £?[r 7/2 ]+£?[p-r 7/2 ]+£?[7r- p] = ^7 + E[p - r l/2 ). (7.37) 



Since 



£[p - r 7/a ] < / £ 

JO 

</ £ 

JO 



L 1 {E fe >o: CTfc < T - 7/2 4 fe Vo} 



fc>0:CT fc <T 7/ / 2 

using moreover Q7.36JI and (|7.37|) . we can estimate 



dt 
dt, 



(7.38) 



q /*oo 

EW t ]dt<\l+ E[ E ^ 

- 70 fc>0:f7 fc <T 7/2 

Since > : < t 7 / 2 }|] = 1 + 7, we obtain 



dt. 



(7.39) 



roo o /-oo 

y o Ew t ]dt < - 7 + (1 + 7 ) y o ^[^ 0) + VW dl 



(7.40) 



Since ^j ^ jumps from n to n — 1 with rate (n + l)n, the expected total time that 



,(0) 



7? 



equals l/((n + l)n), and therefore 
Jo 

Inserting this into (|7.40|) . using (|7.33|) . we arrive at (|7.32|) . 



rt=l 



n=l 



(7.41) 



Lemma 7.11 (The catalyzing function /io,i) One /ias Z/y(/to,i) < ^0,1 / or eac/i 7 > 0. 
Moreover, for each r > 1 and 7 > 0, 



U-y{rh ,i)(x) 
sup — -— < 1. 

*e(o,i] ™ ,i(zJ 



(7.42) 
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Proof Recall that h ,i(x) = h 7 (x) = 1 — (1 — x) 7 (x G [0, 1]). We will show that 

£ (7 ' 0) [^J<7 (7.43) 
for each 7 G (0, 00). The function Xm(l) from Lemma 17.91 satisfies 

Xm(l) = -V-<1 (m>5). (7.44) 
m n 
n=l 

Since Xm{l) 1S concave in 7 and satisfies Xm(0) = 1, it follows that Xm(l) < 1 for all < 7 < 1 
and m > 5. By Lemma 17.101 for all 7 > 1, 

E ^) Woo] < 2 jr\ + \<m (m>7). (7.45) 
fc=i 

Therefore, if m > 7, then m' := E^ m ' ^ [ip'^] < m. It follows by (|7.22j) and Jensen's inequality 
applied to the concave function zwl - (1 — x) z that 

U<yh m (x) < 1 - (1 - x) E{m '° )[ ^ ] = 1 - (1 - x) m ' < h rn {x) (x G [0, 1], 7 > 0). (7.46) 
This shows that h m is £/ 7 -superharmonic for each 7 > 0. By Lemma 17.21 for each r > 1, 



U,(rh m )(x) rU,(h m )(x) 1 - (1 - x) m , , , . 

rM*) " rh m (x) -l-(l-x)- ( 7 - 47 ^ 

By Lemma 17.31 and the monotonicity of , 

U^rh m )(x) ^ 7 (r)(x) l + 7 1 

r/i m (x) -rUi)-l+n(l-(l-i) m ) ^t^ij,. 

Since the right-hand side of Q7.47JI is smaller than 1 for x G (0, 1) and tends to m 1 /m < 1 as 
a; — > 0, since the right-hand side of (|7.48|) is smaller than 1 for x in an open neighborhood of 
1, and since both bounds are continuous, (|7.42|) follows. I 



7.3 Extinction versus unbounded growth 

In this section we show that Lemmas l3.6H3.8l are equivalent to Proposition 13.91 (This follows 
from the equivalence of conditions (i) and (ii) in Lemma 17.121 below.) We moreover prove 
Lemmas 13.61 and 13.81 and prepare for the proof of Lemma 13.71 We start with some general 
facts about log-Laplace operators and branching processes. 

For the next lemma, let E be a separable, locally compact, metrizable space. For n > 0, 
let q n G C+(E) be continuous weight functions, let Q n be continuous cluster mechanisms on E, 
and assume that the associated log-Laplace operators U n defined in ()3.2|) satisfy ()3.3|) . Assume 
that + 1 h G C+(E) is bounded and Z-/ n -superharmonic for all n, let E := {x G E : h(x) > 0}, 
and define generating operators U% : B^ ^(E h ) — ► Biq 11(E) as in ()3.25|) . For each n > 0, 

let (X^ n \x[ n) )he a one-step Poisson cluster branching process with log-Laplace operator U n , 
and let (X^ n) ' h , x[ n) ' h ) be the one-step branching particle system with generating operator U%. 
(In a typical application of this lemma, the operators lA n will be iterates of other log-Laplace 
operators, and Xq , X^ 1 will be the initial and final state, respectively, of a Poisson cluster 
branching process with many time steps.) 
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Lemma 7.12 (Extinction versus unbounded growth) Assume that p G C\^^{E h ) and 
put 

v{x) . = { Kx)p{x) ifxeE h 

P[X) ■ \ ifx£E\E h . [l yj 

Then the following statements are equivalent: 

(i) P*- [\x[ n) > h \ G •] => p(x)5oo + (1 - ?(*))<$„ 

n— *oo 

locally uniformly for x G P ft , 

(ii) i*»[(Af\/i>€-] => e -K^)5 + (i- e -p(^)) 5oo 

locally uniformly for x £ E, 

(iii) U n (Xh)(x) — > p[x) 

n— >oo 

locally uniformly for x G E VA > 0, 

(iv) 30 < Ai < A 2 < oo : U n (Xih)(x) — > p(x) 

n— >oo 

locally uniformly for i£E = 1,2). 

Proof of Lemma 17. 121 It is not hard to see that (i) is equivalent to 

P s *[Thin x (x[ n) ' h ) / 0] — ► p{x) (7.50) 

n— too 

locally uniformly for x G E h , for all < A < 1. It follows from (j3.19j) and l)3.25|) that 
/i(x)P**[Thin A (X 1 n),/l ) / 0] = hU h {\){x) = U(Xh){x) (x G P), so (i) is equivalent to 

(i)' U n (Xh)(x) — p(a:) 

locally uniformly for x £ E V0 < A < 1. 

By (|3.4j) . condition (ii) implies that 

e -U n (Xh)(x) _ E s x r e -\(Xi,h)i y e ~p(x) ( 751 ) 

- 1 n— >oo 



locally uniformly for x G P for all A > 0, and therefore (ii) implies (iii). Obviously (iii) 
(i)' =^(iv) so we are done if we show that (iv)=^(ii). Indeed, (iv) implies that 



E s *[ e - X i( x i n > h ) - e -M x i n — > q (7.52) 



n— >oo 

locally uniformly for x G P, which shows that 

A r„ / / vW 



P fe [c < < C] — ► (7.53) 



for all < c < C < oo. Using (iv) once more we arive at (ii). I 

Our next lemma gives sufficient conditions for the n-th iterates of a single log-Laplace operator 
IA to satisfy the equivalent conditions of Lemma l7.12l Let P (again) be separable, locally com- 
pact, and metrizable. Let q G C+(P) be a weight function, Q a continuous cluster mechanism 
on P, and assume that the associated log-Laplace operator U defined in (|3.2|) satisfies ()3.3|) . 
Let X = (Xq, X\, . . .) be the Poisson-cluster branching process with log-Laplace operator U in 
each step, let ^ h G C+(P) be bounded and W-super harmonic, and let X h = (Xq,X^, . . .) 
denote the branching particle system on E h obtained from X by Poissonization with a IA- 
superharmonic function h, in the sense of Proposition 13.41 
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Lemma 7.13 (Sufficient condition for extinction versus unbounded growth) Assume 
that 

Uh(x) ^ ,_ _ A . 

sup , w < L (7 - 54) 

lee' 1 ' l W 

T/ien i/te process X h started in any initial law C(Xq) G Ai\(E h ) satisfies 

lim |_x£| = 00 or 3k > s.t. x£ = a.s. (7.55) 

fc— >oo 

Moreover, if the function p : E h — > [0,1] defined by 

p( x ) ■= P s *[X%^0 Vn>0] (x G £*) (7.56) 
satisfies ini x&E h p(x) > 0, i/ien /) is continuous. 

Proof of Lemma 17.131 Let A denote the tail event A = {X% / Vn > 0} and let {!Fk)k>o 
be the filtration generated by X h . Then, by the Markov property and continuity of the 
conditional expectation with respect to increasing limits of <7-fields (see Complement 10(b) 
from |Loe631 Section 29] or |Loe781 Section 32]) 

P[X\\ ± Vn > 0|X fc ] = P(A\T k ) — > 1 A a.s. (7.57) 

In particular, this implies that a.s. on the event A one must have P[X^ +1 = 0|X^] — > a.s. 
By (jXTTHl and fl£25D , P 5 -^ / 0] = t/ h l(x) = (Uh(x))/h(x), which is uniformly bounded 
away from one by (ESP. Therefore, = 0\X%] -»■ a.s. on .A is only possible if the 

number of particles tends to infinity. 

The continuity of p can be proved by a straightforward adaptation of the proof of [FS041 
Proposition 5 (d)] to the present setting with discrete time and noncompact space E. An 
essential ingredient in the proof, apart from (|7.54l) , is the fact that the map v ^ P v [X% G • ] 
from M{E) to M\(N(E)) is continuous, which follows from the continuity of Q h . I 

We now turn our attention more specifically to the renormalization branching process X. In 
the remainder of this section, (7fc)fc>o is a sequence of positive constants such that ^2 n j n = 00 
and 7„ — > 7* for some 7* G [0,oo), and X = (X_ n , . . . ,Xq) is the Poisson cluster branching 
process on [0, 1] defined in Section O We put := U ln _ x o • • • o Uy . If / h G C[0, 1] 
is Z// 7fe -superharmonic for all k > 0, then X h and X h denote the branching process and the 
branching particle system on {x G [0, 1] : h(x) > 0} obtained from X by weighting and 
Poissonizing with h in the sense of Propositions 13.31 and 13.41 respectively. 

Proof of Lemma 13.61 Bv induction, it follows from Lemma 17.31 that 

U {n \Xh 1A ) = (A>0 ). (7.58) 

rL= ( i +7fc)-i + j 



00 



It is not hard to see (compare the footnote at (|2.12|l ) that 

OO c 

rj (1 + 7fc) = 00 if an d only if 7^ = 00. (7.59) 

k=0 k=0 

Therefore, since we are assuming that ^ n 7 n = 00, 

W (n) (A^i,i) — ► hi tl , (7.60) 



n— >oo 
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uniformly on [0, 1] for all A > 0. The result now follows from Lemma 17.121 (with h = h\ t i and 

p(x) = i {xe [0,1])). ■ 

Remark 7.14 (Conditions on (7 n )n>o) Our proof of Lemma 13 . 61 does not use that 7„ — > 7* 
for some 7* E [0, 00). On the other hand, the proof shows that ^ n 7„ = 00 is a necessary 
condition for (|3.29J> . O 

We do not know if the assumption that 7„ — ► 7* for some 7* E [0, 00) is needed in Lemma 13.71 
We guess that it can be dropped, but it will greatly simplify proofs to have it around. 

We will show that in order to prove Lemmas 13. 71 and !3.81 it suffices to prove their analogues 
for embedded particle systems in the time-homogeneous processes y 1 (7* E [0, 00)). More 
precisely, we will derive Lemmas 13.71 and 13.81 from the following two results. Below, (U®)t>o 
is the log-Laplace semigroup of the super- Wright-Fisher diffusion y°, defined in 1)3.15(1 . The 
functions £"0,1,7* (7* ^ [O1 00 )) are defined in (|3.34j) . 

Proposition 7.15 (Time-homogeneous embedded particle system with /io,o) 

(a) For any 7* > 0, one has (U^*) n hoQ — ► uniformly on [0,1]. 

' n— >oo 

(b) One has U^hoo — * uniformly on [0, 1]. 

Proposition 7.16 (Time-homogeneous embedded particle system with ^o,i) 

(a) For any 7* > 0, one has (U~*) n (Aho 1) — > pn 1 ~* uniformly on [0, 1], for all A > 0. 

' n^oo ' ' ' 

(b) One has U^(Xhoi) — ► Pn 1 uniformly on [0, 1], for all A > 0. 
Proposition 17.151 fa) will be proved in Section f8. 21 

Proof of Proposition 17.161 (a) By formula (|7.42j) from Lemma 17.111 for each r > 1 the 
function rho,i satisfies condition (|7.54f) from Lemma 17. 131 Set p(x) := P Sx [Yn ,rh °- 1 Vn]. 
Then, by (ETTfll) and 

p(x) = lim P^[Y n 7VVl / 0] = lim ([/" Vl )"l(x) 

n— >oo n—>-oo ' 

- lim ^-r(rh°M > _^feL ^ 
n-»oo rho : i{x) rho,i(x) 

where h\(x) = x (x E [0, 1]) is the Z// 7 *-subharmonic function from Lemma l7.8l It follows that 
m fa;e(o,i] P( x ) > an d therefore, by Lemma 17. 131 p is continuous in x. 

By Lemma 17.131 we see that the Poissonized particle system X^ * 1 exhibits extinction 
versus unbounded growth in the sense of Lemma l7.121 which implies the statement in Propo- 
sition [7^1 (a) . I 

Proof of Propositions 17.151 (b) and 17.161 (b) These statements follow from results in 
FS03|. Indeed, IFS031 Proposition 2] implies that for any / E B + [0, 1] and x E [0, 1], 

U tf^) t -Z^0 if/(0) = /(l) = 0, 

U 0f {x )^Z p * {x ) if/(0) = 0, /(1)>0. (7 - 62) 

To see that the convergence in (|7.62jl is in fact uniform in x E [0, 1] we use the fact that 
each function / E B + [0, 1] with /(0) = /(l) = can be bounded as / < ^l(o,i) f° r some 
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r > 1, and that each function / G -B+[0, 1] with /(0) = and /(l) > can be bounded as 
< / < r l(o,i] f° r some < e < 1 and r > 1. Therefore, by the monotonity of , it 
suffices to show that £/ t °(rl( 0i i)), W t (rT(o,i]), and Z^ t °(el{i}) converge uniformly on [0,1]. By 
FS03, Lemma 15], these functions are continuous for each t > 0, and since moreover the limit 
functions are continuous, it suffices to show that the convergence is monotone. Thus, we claim 
that 

W°(rl (0)1) )i0 (r>l), 

z4V1(o,i]Hp5,i,7* (*■>!)> ( 7 - 63 ) 

W°( e l{i})TPo,i, 7 * (0<e<l). 
By (an obvious analogue of) Lemma I7.2( it suffices to show that l(o,i) an d l(o,i] are M?~ 
super harmonic, while l^j is Z^ t °-subharmonic for each t > 0. Let (y/ 1 ' 1 )t>o be the branching 
particle system obtained from (3^)t>o by Poissonization with the constant function h\ \ := 1. 
Then y ^ 1 * 1 is a system of binary splitting Wright-Fisher diffusions, which was also studied 
in [ESSE]. One has (compare (TTTTTl) ^ 

l4l m (x) = P[Thm u o 1(o i) (5 x ) * 0] = P^[Thin 1(01) (y t °' ft1 ' 1 ) + 0] = P s *[Y t ° M >\(0, 1)) > 0]. 

(7.64) 

Likewise, 

^°l (0il] (x) = P^[y t °^ 1 ((0,l])>0] and U t °l {1} (x) = P s *[Y t ° M >\{l})>0}. (7.65) 

Using the fact that the points 0, 1 are traps for the Wright-Fisher diffusion and that in a binary 
splitting Wright-Fisher diffusion, particles never die, it is easy to see that P s *[Y t ' M ((0, 1)) > 
0] and P s *[Y t °' hl - 1 ((0, 1]) > 0] are nonincreasing in t, while P s * [Y t °' hlA ({1}) > 0] is nonde- 
cr easing in t. I 

We now show that Propositions 17.151 and ITTol imply Lemmas 13.71 and 13.81 respectively. 

Proof of Lemma 13.71 We start with the proof that the embedded particle system X h °'° is 
critical. For any / G B + [0, 1] and k > 1, we have, by Poissonization (Proposition 13.41) and the 
definition of X, 

h Qfi {x)E~ k ^[(x h _Znf)] = £- fe >APois(Vo<y) [( ^M> = £r fe ^[<Pois(v*-fc+i),/>] 

= E- k > 5 *[{X-k + iMM = + l)E[(Z2,h , f)] = (i + i)<r?-W/>, 

(7.66) 

where is the invariant law of y2 from Corollary 15.41 In particular, setting / = 1 gives 

/i ,o(*)£- fcA [|^° fc +il] = M*) b y <E2SD- 

To prove (|3.3Uj) . by Lemma 17.121 it suffices to show that 

W(*°(Afco,o) — > (7.67) 

uniformly on [0, 1] for all < A < 1. We first treat the case 7* > 0. Then, by Theorem 13 . 21 (a) . 
for each fixed I > 1 and / G C+[0, 1], 

Uyn-iO-'-oUrn-if ^(VrYf ( 7 - 68 ) 
uniformly on [0, 1]. Therefore, by a diagonal argument, we can find l{n) — > 00 such that 

IKz4*y (n) Vo -^_x ° • • • °Wy tt _, (n) VolL _-± °- ( 7 - 69 ) 
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Using the fact that the function /io,o is £/ 7 -superharmonic for each 7 > and the monotonicity 
of the operators we derive from Proposition 17.151 fa) that 

U^(Xh ,o) < U Yn _ 1 o • • • oW 7 „ j/io.0 — (7.70) 



n— >oo 



uniformly on [0, 1] for all < A < 1. This proves (j7.67j) in the case 7* > 0. 

The proof in the case 7* = is similar. In this case, by Theorem 13.21 (b), for each fixed 
t > and / G C+[0,1], 

Mrn-a • • -"^w/W — *470*0 Vx„ -> x G [0, 1], (7.71) 

which shows that W 7n _ 1 ° • • • o W 7fc , t) / converges to f uniformly on [0, 1]. By a diagonal 
argument, we can find t(n) — > 00 such that 

\\U°(h , ) -Uy^ o • • • oZ^ (t(n)) (ML n ~=^ °> ( 7 - 72 ) 

and the proof proceeds in the same way as before. I 

Proof of Lemma 13.81 By Lemma 17.121 and the monotonicity of the operators it suffices 
to show that 

(i) limsupW< n >(fc ,i)<Po.i.«v J 

(7.73) 



(i) limsupW( n )(V)<Po,i, 7 *, 

(ii) limmfW( n )(±Vi)>^ 1)7 *, 



uniformly on [0,1]. We first consider the case 7* > 0. By (j7.68f) and a diagonal argument, we 
can find l(n) — ► 00 such that 

IKw 7 *y ( "V -Uf^.i o ■ ■ ■ ow^^oaL ^ °- ( 7 - 74 ) 

Therefore, by Proposition 17.161 (a), the fact that /io,i is lL k -super harmonic for each k > 0, 
and the monotonicity of the operators IL, we find that 

U [n) h 0) i < W in _ x o ■ ■ -oUy n _ Kn) hQ^ — + Po,i, 7 *> (7-75) 



uniformly on [0, 1]. This proves (|7.73|) (i)- To prove also (|7.73|) (ii) we use the W 7 -su&harmonic 
(for each 7 > 0) function h\ from Lemma l7.8l By Lemma l7.2l also hh\ is ^ 7 -subharmonic. By 
bounding from above and below with multiples of /io,i it is easy to derive from Proposi- 
tion [7^1 (a) that 

[UfT&hi) -— pS )1>r (7.76) 
uniformly on [0, 1]. Arguing as before, we can find l(n) — > oo such that 

||(^.) J(n) (^-^-i°"- o ^n-Kn)(3MlL f ^ ) 0. (7.77) 

Therefore, by (|7.76f) and the facts that \h\ is £/ 7fc -subharmonic for each k > and ^/ii < |/io,i) 

uniformly on [0, 1], which proves (|7.73|) (ii). The proof of (|7.73|) in case 7* = is completely 
analogous. I 
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8 Extinction on the interior 

8.1 Basic facts 

In this section we prove Proposition 17.151 (a). To simplify notation, throughout this section 
h denotes the function /io,o- We fix < 7* < 00, we let Y h := Y" 1 > h denote the branching 
particle system on (0, 1) obtained from y" 1 = (y^ , yj , . . .) by Poissonization with h in the 
sense of Proposition 13. 4( and we denote its log-Laplace operator by Ub,. We will prove that 

p{x) := P S * [y£ / Vn > 0] = (s € (0, 1)). (8.1) 

Since for each n fixed, x 1— ► p n (x) := P &X \Y^ ^ 0] is a continuous function that decreases to 
p(x), <|8.1|) implies that p n (x) — > locally uniformly on (0, 1), which, by an obvious analogon 
of Lemma IY. 121 yields Proposition 17. 151 fa). 
As a first step, we prove: 

Lemma 8.1 (Continuous survival probability) One has either p{x) = for all x G (0, 1) 

or there exists a continuous function p : (0, 1) — > [0, 1] such that p(x) > p{x) > for all 
x 6 (0,1). 

Proof Put p(x) := h(x)p(x). We will show that either p = on (0,1) or there exists a 
continuous function p : (0, 1) — ► (0, 1] such that p > p on (0, 1). Indeed, 

p(x) = h{x)P 6x \Y% / Vn > 0l = lim h(x)P 5 * \Y% / 0l 

n—>oo 1 f8 2) 

= h(x) lim (Ub) n l(x) = lim (U^) n h(x) (z € (0,1)), 

n— »oo ' n— >oo 

where we have used (|3,19[) and H3.25() in the last two steps. Using the continuity of IA^* with 
respect to decreasing sequences, it follows that 

U r p = p. (8.3) 

We claim that for any / G £?r u[0, 1], one has the bounds 

(T2, f) < Urf(x) < (1 + 7)(r^ /) (7 > 0, x G [0, 1]). (8.4) 

Indeed, by Lemma 1731 U 7 f(x) > 1 — E[(l — /(y2(0)))] = (TZ,f), while the upper bound in 
(H3J) follows from TTJU . 

By Remark 15.51 (0,1) 9 x 1— > (rj, /) is continuous for all / G S[ 0)1 ][0, 1]. Moreover, 
(T2,f) = for some x G (0,1) if and only if / = almost everywhere with respect to 
Lebesgue measure. 

Applying these facts to / = p and 7 = 7*, using (|8.3jl . we see that there are two possibilities. 
Either p = a.s. with respect to Lebesgue measure, and in this case p = by the upper bound 
in (j8,4[) . or p is not almost everywhere zero with respect to Lebesgue measure, and in this case 
the function x 1— > p(x) := (rj, /) is continuous, positive on (0, 1), and estimates p from below 
by the lower bound in (|8.4|) , I 
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8.2 A representation for the Campbell law 

(Local) extinction properties of critical branching processes are usually studied using Palm 
laws. Our proof of formula (|8.1|) is no exception, except that we will use the closely related 
Campbell laws. Loosely speaking, Palm laws describe a population that is size-biased at a 
given position, plus 'typical' particle sampled from that position, while Campbell laws describe 
a population that is size-biased as a whole, plus a 'typical' particle sampled from a random 
position. 

Let V be a probability law on AA(0, 1) with Jj^/ Q j\ 7'(d^)|f| = 1. Then the size-biased law 
Vsize associated with V is the probability law on AA(0, 1) defined by 



?W-):= / V{du)\v\lr y (8.5) 

The Campbell law associated with V is the probability law on (0, 1) x AA(0, 1) defined by 

P C amp(A x B) := ! V(du) u(A)lr B} (8.6) 

for all Borel-measurable A C (0,1) and B C Af(0, 1). If (v, V) is a (0,1) x AA(0, l)-valued 
random variable with law "Pcamp, then CiV) = V s i zc , and v is the position of a 'typical' 
particle chosen from V. 
Let 

P x ' n (-):=P s *[Yj?e-] (8.7) 

denote the law of Y h at time n, started at time with one particle at position x G (0, 1). 
Note that by criticality, J^ ^ 'P :c ' n (di/)|^| = 1. Using again criticality, it is easy to see that 
in order to prove the extinction formula ()8.1|) . it suffices to show that 

Jim P*£({1,...,N})=0 (*G(0,1), N>1). (8.8) 

In order to prove ()8.8|) . we will write down an expression for VQ^ mp - Let Q h denote the offspring 
mechanism of Y , and, for fixed x G (0, 1), let <3camp( x > ' ) denote the Campbell law associated 
with Q h (x, • ). The next proposition is a time-inhomogeneous version of Kallenberg's famous 
backward tree technique; see |Lie81| Satz 8.2]. 

Proposition 8.2 (Representation of Campbell law) Let (vfc,T4)fc> be the Markov pro- 
cess in (0, 1) x M(0, 1) with transition laws 

P[(v fe+ i, V k+1 ) G • | (v fc , V k ) = (x, u)} = Q h Camp (x, • ) ((x, v) G (0, 1) x Af(0, 1)), (8.9) 

started in (vo,Vo) = (5 X , 0). Let (y ft >( fc )) fc ^ 1 5e branching particle systems with offspring 
mechanism Q h , conditionally independent given (v/%, T4)fc>o, started in Yn = Vk — 5 Vk . 
Then 

n 

-P X cL v = c{^5 Vn+ Y. Y n^)- (8-10) 

k=l 

Formula (|8.1()|) says that the Campbell law at time n arises in such a way, that an 'immortal' 
particle at positions vo, . . . , v n sheds off offspring V\ — <5 Vl , . . . , V n — 5 Vn , distributed according 
to the size-biased law with one 'typical' particle taken out, and this offspring then evolve 
under the usual forward dynamics till time n. Note that the position of the immortal particle 
( v k)k>o is an autonomous Markov chain. 

We need a bit of explicit control on Qc amp - 
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Lemma 8.3 (Campbell law) One has 

Q h Camp (x,Ax B) = 2±1 J P[Poi S (hz£) e d X ]x(.A)l {xeA} , (8.11) 



where the random measures Z2 are defined in \3. 

Proof By the definition of the Campbell law (|8.6|) . and l|3.24|) . 

Q^ mp (x,Ax B) = J Q h (x,d X )x(A)l {xeB } 



h{x) 

Recall that by (TTTl) . 



J P[Pois(^ ) G dx]x(^)l {xeB} + (l - ^7^) • 0. 



o 



6 yT{ _ t/2) dt, (8.13) 



where (y2 (i))teR is a stationary solution to the SDE (|3.6|) with 7 = 7*. By Lemma EH3 the 
transition law of the Markov chain (yk)k>o from Proposition 18.21 is given by 

P[v fc+1 G dy|v fc = x] = ^^£[Pois(fcZg*)(dy)] = ±+£-h(y)Tf (dy), (8.14) 

where is the invariant law of y2 from Corollary 15.41 In the next section we will prove the 
following lemma. 

Lemma 8.4 (Immortal particle stays in interior) The Markov chain (v/J^x) started in 
any vq = x £ (0, 1) satisfies 

( v fc)fc>o has a cluster point in (0,1) a.s. (8.15) 

We now show that Lemma f8.4| together with our previous results, implies Proposition l7. 151 (a). 
Proof of Proposition I7T15I (a) We need to prove (|8.1[) . By our previous analysis, it suffices 



to prove (|8.8|) under the assumption that p ^ 0. By Proposition 18.21 

n 



k=l 



Conditioned on (vfc,14)fc>o, the (Y \ )k=l,...,n are independent random variables with 

p[Y^f + 0] > P[Y*M + Vm > 0] = P[Thin,(14 - S v J + 0]. (8.17) 

Therefore, (|8,8|) will follow by Borel-Cantelli provided that we can show that 

00 

^P[Thin p (V fe -<5 v J^0| Vfc-i] = 00 a.s. (8.18) 

fe=l 
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Define f(x) := P[Thin p (V fc - 5 Vk ) / 0|v fe _i = x] (x € (0,1)). We need to show that 
^2h=i f( x ) = 00 a - s - Using Lemma l8~TI and Lemma Is73l we can estimate 

f(x) > P[Thmp(V k -5 Vk )^ 0|v fe _! = x] = [ Q h Climp (x,dy,du){l - (1 - p)^} > 

(8.19) 

for all x £ (0,1). Since Q 7 *, defined in ([3.8 j) . is a continuous cluster mechanism, also 
( 3camp( x ' ') * s continuous as a function of x, hence the bound in (|8.19|) is locally uniform 
on (0, 1), hence Lemma 18 . 41 implies that there is an e > such that 

P[Thm p (V k - <5 V J + 0|v fc _i] > e (8.20) 

at infinitely many times k — 1, which in turn implies ()8.18j) . I 

8.3 The immortal particle 

Proof of Lemma 18.41 Let K(x,dy) denote the transition kernel (on (0,1)) of the Markov 
chain (v fc ) fc > , i.e., by (|8.14|). 



K(x,dy) = (1 +7 *)MLJ/W(dy). (8.21) 
x(l — x) 



It follows from lfB~2^)) that 



v > i/m ^ + 2 7 *)(i + 37*) v ; 



Set 



Then 



g(x) := J K(x, dy)y(l - y) - x(l - x) (x € (0, 1)). 



5.23) 



n-l 



M n := v n (l - v n ) - ^ 5 (v fc ) (n > 0) (8.24) 

k=0 

defines a martingale (M n ) n >o- Since g > in an open neighborhood of {0, 1}, 

-P[( v fc)fc>o has no cluster point in (0, 1)] < P\ lim M n = —00] = 0, (8.25) 

— n^oo 

where in the last equality we have used that (M n ) ra >o is a martingale. I 



9 Proof of the main result 

Proof of Theorem II . 41 Part (a) has been proved in Section 1531 It follows from (|2.12|) . (|2.13() . 
(|.3.10[) . and (f3. lift that part (b) is equivalent to the following statement. Assuming that 



00 

(1 



Eln = 00 and (ii) 7„ — > 7* (9.1) 
n. — »no 



n— >oo 
71=1 
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for some 7* G [0, 00), one has, uniformly on [0, 1], 

W 7n _ lO ...o^ (p) -^p* r7 ,, (9.2) 

where p ; * ~* is the unique solution in Hi r of 

(i) U^p* =p* if < 7* < 00, 

(ii) \ x (l- x )£ I p*( x )-p*{ x ){l-p*( x )) = Q (x€[0,l]) if 7 * = 0. (9 ' 3) 

It follows from Proposition 13,51 that the left-hand side of (|9.2|) converges uniformly to a limit 
P* ri * which is given by ()3.34j) . We must show 1° that p* r7 * £ Hi jT and 2° that p1 ri * is the 
unique solution in this class to (j9.3|h We first treat the case 7* > 0. 

1° Since Po,o,7* = and p* l * = 1, it is obvious that ^0,0,7* e ^0,0 and p* i ~* £ 
Therefore, by symmetry, it suffices to show that Pq i~* £ Ho,i- By Lemmas 17.81 and 17.111 
x < p < 1 — (1 — x) 7 implies x < U lk p < 1 — (1 — x) 7 for each k. Iterating this relation, using 
(|9~2l . we find that 

x<p* oxr (x) < 1- (1-x) 7 . (9.4) 

By Proposition I5.ll! the left-hand side of (|9.2[) is nondecreasing and concave in x if p is, so 
taking the limit we find that p$ 1 7 * is nondecreasing and concave. Combining this with Q9.4|) 
we conclude that Po,i,7* ^ s Lipschitz continuous. Moreover p^i^. (0) = and Pq 17 *(1) = 1 so 
^0,1,7* G ^0,1 • 

2° Taking the limit n — ► 00 in (U^*) n p = U~ / *(U~ / *) n ~ 1 p, using the continuity of £/ 7 * 
(Corollary EHOJ) and (|9~2)) . we find that W 7 *p* « = p* r7 «- It follows from (|9~2j) that p^* r » is 
the only solution in 7i; r to this equation. 

For 7* = 0, it has been shown in |FS031 Proposition 3] that p ; * r0 is the unique solution 
in Hi r to ()9.3|) (ii). In particular, it has been shown there that Pq 1 ^ s twice continuously 
differentiable on [0, 1] (including the boundary). This proves parts (b) and (c) of the theorem.! 



A Appendix: Infinite systems of linearly interacting diffusions 
A.l Hierarchically interacting diffusions 

For any N > 2, the hierarchical group with freedom N is the set f2jv of all sequences £ = 
£2; • • • )> w ith coordinates in the finite set {0, . . . , N — 1}, which are different from only 
finitely often, equipped with componentwise addition modulo N. Setting 

||f|| := min{n > : £ k = VA: > n} (£ G Sl N ), (A.l) 

||£ — r; 1 1 is said to be the hierarchical distance between two sites £ and r/ in Q]\f. 

Let D C M. d be open and convex, and let W be a renormalization class on D. Let cr 
be a continuous root of a diffusion matrix w G W as in Remark 11.21 Consider a collection 
x = (x^)^ e fi N of D-valued processes, solving a system of SDE's of the form 

00 

dx ?W = E / ^(4 +1 W- x ?W) dt + ^ CJ ( x €W) dB €W (^°> ( A - 2 ) 

fc=0 
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where (B^)^q n is a collection of independent standard Brownian motions, with initial condi- 
tion 

x ( (0)=fl£D (^en N ). (A.3) 

Here the (ck)k>o are positive constants satisfying Yl,k c k/N k < oo, and x|(i) denotes the 
k-block average around £: 

(Note that |{ry : ||£ — r/|| < fc}| = N k .) Under suitable additional assumptions on a, one can 
show that (|A.2|) has a unique (weak or strong) solution (sec DG93a | IDG961 ISwaOOj ) . We call 
x a system of hierarchically interacting D-valued diffusions with migration constants (ck)k>o 
and local diffusion rate Wij = X^fc °~ikO~jk- Such systems are used to model gene frequencies or 
population sizes in population biology |SF83j . 

The long-time behavior of the system in (|A.2|) depends crucially on the recurrence versus 
transience of the continuous-time random walk on fijy which jumps from a point £ to a point 
7] 7^ £ with rate 

oo 

*(v -*)--= E w&- ( A - 5 ) 



N 2k- 



This random walk is recurrent if and only if 



oo 



— = oo, where dh '■= fc+n (A. 6) 

d k ^ N n v ' 

k=0 K n=0 



(see IDG93al lKle96 ; a similar problem is treated in DE68 ). Assuming that the law of x(i) 
converges weakly as t — > oo to the law of some D ^-valued random variable x(oo), one expects 
that in the recurrent case x(oo) must have the following properties: 

(i) x € (oo) = x^(oo) a.s. V^r/Gfijv, 

(ii) x $ (oo) E a w D a.s. V££Q N . 1 j 

Here d w D is the effective boundary of D, defined in (|2.3j) . If x(i) converges in law to a limit 
x(oo) satisfying (|A.7|) . then we say that x clusters. In the transient case, it is believed that 
solutions of (|A.2|) do not cluster. (For compact D these facts were proved in |SwaOO| .) 

An important tool in the study of solutions to (|A.2j) is the so-called interaction chain. 
This is the chain (x[J(i), Xq(£), . . .) of block-averages around the origin. Heuristic arguments 
suggest that in the local mean field limit N — > oo, the interaction chain converges to a certain 
well-defined Markov chain. 

Conjecture A.l Fix w S W, 9 £ D, and positive numbers (ck)k>o such that for N large 
enough, ^2k c k/N k < oo. For all N large enough, let x. N be a solution to L4.jjj) - L4, 3\) . and 
assume that tN are constants such that, for some n > 1, limAr^oo N~ n tN = T € [0, oo). Then 

(4> n (t N ),...,4>°(t N )) =► (/-„,..., jy), (a.s) 

\ / N— >oo 

where (PH n , ... is a Markov chain with transition laws 

P[r k € dy|/Vi =*] = ^ F(k)w (dy) (x G A < fc < n - 1) (A.9) 
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and initial state 

r n = YT, where dy t = c n (9 - y t )dt + V2^ n \y t )dB t , y = 9, (A.10) 
and is a root of the diffusion matrix F^w. 

Rigorous versions of coni ect ure I A . 1 1 have been proved for renormalization classes on D = [0, 1] 
and I) = [0,oo) in jl )C93al IT)093hj . 

Note that the iterated kernels K w, ^ n ' defined in Q2.4j) are the transition probabilities from 
time — n to time of the interaction chain in the mean-field limit: 

K^ n \dy) = P[I$ G dy\I™ n = x] (xeD, n> 0). (A.ll) 

Lemma [2 . 6l expresses the fact that the system clusters in the local mean-field limit N — > oo. 
The condition s n —* oo in Lemma |2.6I means that X}fc>o c = °°> wmc h> in a sense, is the 
N —* oo limit of condition (IA.6I). 



A. 2 The clustering distribution of linearly interacting diffusions 

Let -D C M d be open, bounded, and convex, and let W be a renormalization class on D. 
Fix migration constants (cfc)fc>o and assume that s n — ► oo and s n+ i/s n — > 1 + 7* for some 
7* € [0,oo]. Recall the definition of the iterated probability kernels K w ^ n ' in (|2.4|) . Recall 
Conjecture 12. 71 Assuming that the rescaled renormalized diffusion matrices s n F^w converge 
to a limit it;*, we can make a guess about the limit of the iterated probability kernels K w ^ n \ 

Conjecture A. 2 (Limits of iterated probability kernels) Assume that s n F^ n 'w — > w* 
as n —* 00. Then, for any uu £ W, 

n— >oo 

where K* has the following description: 
(i) 7/ < 7* < 00, i/ien 

K* = lim P x [P n ' G •], (A.13) 

n— >oo 

where (In )n>o Markov chain with transition law P[I^ l+ i G ■ \In = x] = v 1 /" 1 > w . 



(ii) If 7 * = 0, then 



K* = lim P*[J° G •], (A.14) 

1^00 



where (J°) s >o is i/ie diffusion process with generator ]T\ j=\ w*j (y) Qy.Q y . ■ 
(iii) 7/7* = oo, £/ien 

= lim i^ 7 '™*. (A. 15) 

7^00 

For each N > 2, let x w = (x^ ^gp^ be a system of hierarchically interacting diffusions as in 
(|A.2|) and (|A.3|) . If 7* = 0, then because of Conjectures IA.1I and IA.21 we expect 3 that 

lim lim £(x^(A n T)) = K* e (T > 0), (A.16) 

71— »oo N^oo 



3 For 7* > 0, the situation is more complex. In this case at the right-hand side of l|A.160 we expect the law 
fjjP [yT S <lx]K*, where y solves the SDE dy t = -^-(8 — Vt)dt + y^a" (y t )dB t and a* is a root of the diffusion 
matrix w* . Note that in this case the right-hand side of HA.16I depends on T. 
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where K* is the kernel in l|A.14|) . 

In particular, consider the case that the migration constants (ck)k>o are of the form Ck = r k 
for some r > 0. In this case, s n+ i/s n — ► £ V 1, and s n — > oo if and only if r < 1. One can 
check (see l|A.6|l ) that for fixed N > 2, the random walk with the kernel a in (|A.5[) is recurrent 
if and only if r < 1. The critical case r = 1 corresponds to a critically recurrent random walk. 
For a precise definition of critical recurrence, see |Kle 96 , formula (1.15)]. For r = 1, we expect 
that the double limit in (jA.16j) can be replaced by a single limit. More precisely, for each fixed 
N > 2, we expect that 

lira C(^(t)) = Kg. (A.17) 

t— >oo 

In this case, we call Kg the clustering distribution of x JV . The clustering distribution of linearly 
interacting isotropic diffusions was studied in |Swa00| . We expect (|A.17|) to hold, even more 
generally, for all systems of linearly interacting diffusions with a critically recurrent migration 
mechanism. In particular, we expect (|A,17|) to hold for symmetric nearest-neighbor interaction 
on Z d in the critical dimension d = 2. If one is ready to make this enormous leap of faith, 
then combining Conjectures 12.71 and lA.2l one arrives at the following conjecture. 

Conjecture A. 3 (Critical clustering) Let D C M. d be open, bounded, and convex, and let 
W be a renormalization class on D. Assume that the asymptotic fixed point equation \2.1b)) (ii) 
has a unique solution w* in W. Let a be a continuous root of a diffusion matrix w G W. Let 
x = (x^)^ eZ 2 be a D -valued process, solving the system of SDE's 

dxf(t)= Y, (x ?? (t)-x 5 (t))dt + a(x € (t))d J B € (t), (A.18) 

with initial condition x^(0) = 9 G D (£ G Z 2 ). Then 

Xf(t)=>i«, (£ez 2 ), (A.19) 

t^oo 

where (/f ) s >o * s the diffusion with generator Yli j w ij(v) dy^dy an( ^ ^ n ^ a ^ condition I® = 9. 
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